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ABSTRACT 

Steam condensation heat transfer measurements were made 
in a 5-tube test condenser having an additional perforated 
tube to simulate up to 30 active tubes. Results were 
obtained for smooth tubes and roped tubes wrapped with wire. 
A Sieder-Tate equation was used to correlate the inside 
heat-transfer coefficient. For smooth tubes, a leading 
coefficient of 0.029 was found, while it was 0.061 for the 
roped tubes. The average condensing coefficient measured 
for 30 smooth tubes was 0.59 times the Nusselt coefficient 
calculated for the first tube. When the smooth tubes were 
wrapped with wire, this ratio increased up to 0.86. 

Further, roped tubes without wire experienced a ratio of 
0.63, while roped tubes wrapped with wire resulted in a 
ratio of 0.86. These preliminary data show that wire- 
wrapped tubes may lead to a significant reduction in 
condenser surface area. 
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NOMENCLATURE 



2 

A Outside, heat-transfer area of one tube (m ) 

o ' 

2 

A^ Inside, heat-transfer area of one tube (m ) 
Sieder-Tate coefficient 

Specific heat of water evaluated at Tj^ (KJ/kg*K) 

0 14 

Ci Correction factor (u/u ) 
f ^ 

Inner diameter of the tube (m) 

D Outer diameter of the tube (m) 

o 

2 

g Acceleration of gravity (9.81 m/s ) 

h. Experimentally-determined value for the inside, heat- 
^ transfer coefficient (W/m^’K) 

hgg Latent heat of vaporization (KJ/kg) 

hj^ Local, outside, heat-transfer coefficient for the Nth 

^ tube (W/m^K) 

h Heat-transfer coefficient calculated from the Nusselt 

equation (W/m^-K) 

h. Outside heat-transfer coefficient for the first tube 
( W/m2 • K ) 

k^ Thermal conductivity of the condensate film (W/m*K) 

k Thermal condictivity of the cooling water evaluated 

at Tj^ (W/m*K) 

k^ Thermal conductivity of titanium (W/m*K) 

L Condensing length (m) 

LMTD Logarithmic Mean Temperature Difference (°C) 
m Slope of the least-squares-f it , straight line 

m Mass flow rate of cooling water (kg/min) 
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N The munber of tubes in a column or the tube number of 

a given tube 

Nu Water-side Nusselt number 

Prandtl number evaluated at Tj^' 

Q Heat transfer rate (W) 

2 

q” Heat flux based on outside area (W/m ) 

Rg Water-side Reynolds number 

2 

R^ Fouling thermal resistance (m K/W) 

R Wall thermal resistance based on the outside area 

^ (m^K/W) 

R^ Outside, local, heat-tranfer coefficient ratio (h^^/hj^) 
R^ Outside, average, heat-tranfer coefficient ratio 

(hj^/hi> 

S/D Spacing-to-diameter ratio of tubes 

Average cooling water bulk temperature (°C) 

Cooling water inlet temperature (°C) 

Cooling water outlet temperature (°C) 

Wall temperature (°C) 

Average condensate film temperature (°C) 

T . Saturation temperature of steam (°C) 

s ^ tz 

Vapor (steam) temperature (°C) 

2 

Overall heat-transfer coefficient (m K/W) 

2 

Outside heat-transfer coefficient (m K/W) 

Cooling water velocity (m/s) 

X Sieder-Tate parameter (X = Re Pr- ) (— ) 

^w 

AT Temperature difference (°C) 
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Greek Symbols 



? Heat capacity parameter (C AT/h^ ) 
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H Dynamic viscosity of water evaluated at Tj^ (N*s/m ) 
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Dynamic viscosity of condensate evaluated at (N*s/m ) 
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u Dynamic viscosity of water evaluated at T (N*s/m ) 

C Acceleration parameter (KAT/ph^g) 
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I. HISTORICAL BACKGROUND 



In recent years, there has been a continued interest in 
the reduction of the size and weight of propulsion systems 
aboard both surface vessels and submarines. Especially, 
actual dimensions of naval condensers have a critical bear- 
ing on cost and performance of the ship. Often, compactness 
is more important than thermal effectiveness when the over- 
all performance of the ship is considered; in a submarine, 
the diameter of the pressure hull can depend on the dimen- 
sions of the condenser. 

The importance of compactness justifies measures to 
raise the overall heat-transfer coefficient of condenser 
tubes despite the penalties which may occur, i.e., the 
increased pumping power and tube cost. 

Naval condenser design is based upon the Heat Exchange 
Institute (HEI) specifications for steam condensers [Ref. 1] 
and also the standards of the Tubular Exchange Manufacturers 
Association (TEMA) [Ref. 2]. Search [Ref. 3] investigated 
the present condenser design processes, including the 
feasibility of enhanced heat transfer in naval condensers. 

He concluded that the current design is very conservative, 
and he predicted that a forty-percent reduction in condenser 
weight and volume could be achieved depending on the heat- 
transfer enhancement method used. 
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In recent years, many research efforts have been 
directed to the study of heat-transfer enhancement tech- 
niques and their application to heat-exchanger design. Webb 
[Ref. 4] has summarized extensive works of augmentation 
techniques. At the Naval Postgraduate School, Beck [Ref. 

5], Pence [Ref. 6], Reilly [Ref. 7], Fenner [Ref. 8] and 
Ciftci [Ref. 9] conducted experimental research into various 
kinds of enhancement schemes employing a single-tube test 
condenser. 

The above-mentioned investigations concluded that, for 
the same diameter tube, the overall heat-transfer coeffi- 
cient of enhanced tubes can exceed those for smooth tubes by 
almost 100 percent. Reilly and Fenner [Ref. 10] revealed 
that most of the above mentioned augmentation occurred on 
the cooling-water side due to a combination of increased 
surface area, and increased turbulence and swirl in the 
cooling water flow. Little or no improvement occurred on 
the steam side. Eissenberg [Ref. 11] performed an extensive 
study on condenser-tube , heat-transfer coefficients using a 
multi-tube bundle. 

In order to investigate the outside heat-transfer 
performance of various enhanced tubes in tube bundles, 
research was conducted at the Naval Postgraduate School. 
Noftz [Ref. 12] modified a test apparatus initially designed 
by Morrison [Ref. 13] to simulate a tube bundle using five 
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active tubes arranged in a vertical plane. A perforated 
tube was located at the top of the bundle, through which 
water was flooded to simulate bundles having up to 30 tubes 
in a vertical row. His investigations determined that the 
heat-transfer coefficients for a given tube of the tube 
bundle increased as the mean vapor velocity increased, but 
decreased as the amount of condensate inundation increased. 
The experimentally found values for the heat-transfer 
coefficients were comparative with the Nusselt theory. 

Based upon research done currently, it is evident that 
present day smooth-tube steam condensers, operating under 
typical conditions, have limitations in their thermal 
efficiency, due to a large thermal resistance which occurs 
on the tube side of the condenser. This resistance is 
generally larger than any of the thermal resistances that 
occur on the steam side, in the tube wall, those due to 
fouling, or due to noncondensable gases. However, employing 
enhanced tubes, the inside heat-transfer coefficient can be 
increased by 100 to 200 percent over the smooth-tube case. 
The outside heat-transfer coefficient, on the other hand, is 
increased by only 10 to 50 percent. In this situation, the 
thermal resistances on the inside and outside of the tube 
can be approximately equal. 
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Webb [Ref, 4] reported that the dominant thermal resis- 
tance in film condensation is that of conduction across the 
condensate film and, therefore, a surface geometry that 
promotes reduced film thickness will provide enhancement. 

Thomas, et al, [Ref, 14] tested ammonia condensation on 
a smooth tube with a wire wrapped in a helical manner. The 
measured condensing coefficient was approximately three 
times that predicted by the Nusselt equation for a smooth 
tube. Surface-tension forces draw the condensate to the 
base of the wires, which act as condensate run off channels. 

Webb [Ref. 4] stated that, when noncondensables are 
present, an additional thermal resistance is introduced in 
the gas at the vapor-liquid interface. Mixing in the gas 
film will substantially reduce this thermal resistance. 
Therefore, the maintenance of high vapor velocities, or 
special surface geometries that promote a higher heat- 
transfer coefficient in the gas film will substantially 
alleviate the performance deterioration due to 
noncondensables . 

Cunningham [Ref. 15] presented in his paper that, for 
the roped tubes on the vapor (shell) side, the enhancement 
is achieved by improved condensate drainage, while on the 
coolant (tube) side the helical ridges increased turbulence 
and, as a result, the inside convective coefficient. 
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Improvements on the condensing side up to 100 percent 
have been reported for single-tube tests [Refs. 16;17] . 
Although titanium has a low thermal conductivity, it pro- 
vides a high resistance to erosion and water-side fouling. 
Titanium tubes with enhancement both inside and outside are 
commercially available through the Wolverine Tube Division 
of Universal Oil Products, Inc. The applications of these 
tubes having all the inherited properties of titanium are 
promising for naval condensers. 

The goals of this thesis were therefore to; 

1. Obtain baseline heat-transfer performance data for 
the test condenser utilizing 16 mm O.D. smooth 
titanium tubes. 

2. Conduct steam condensation tests with the following 
enhanced tube geometries to determine steam-side 
heat-transfer coefficients in relation to smooth- 
tube performance; 

a. Wolverine "roped" tubes. 

b. Wolverine "roped" wrapped with titanium wire. 

c. Smooth tubes wrapped with titanium wire. 



19 



II. THEORETICAL BACKGROUND 



The combined effect of vapor shear and inundation on the 
condensate film heat-transfer coefficient for cylindrical, 
horizontal tubes within tube bundles is a very complex and 
still insufficiently-understood subject, which is of impor- 
tance to the efficient design of steam condensers. Although 
many researchers have studied this subject both theoretic- 
ally and experimentally, there is no accurate methodology 
available for predicting the condensate-film, heat-transfer 
coefficient within tube bundles. 

In 1916, Nusselt conducted his pioneering analysis for 
the simple case of condensation occurring on the outside of 
a single, isolated, horizontal tube. He idealized the 
problem by making the following assumptions for single tubes 
as stated by Nobbs [Ref. 18]. 

1. The wall temperature is constant. 

2. The flow is laminar in the condensate film. 

3. Heat transfer in the condensate is by conduction, 
and subcooling may be neglected. 

4. The fluid properties are constant within the 
condensate film. 

5. The forces due to hydrostatic pressure, surface 
tension, inertia, and vapor-liquid interfacial 
shear are negligible when compared to the viscous 
and gravitational forces. 

6. The surrounding steam and vapor/liquid interface 
are at saturation temperature. 
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7. The film thickness is small when compared with 

normal tube diameters and the effects of curvature 
are small. 



Based on the above assumptions, Nusselt predicted the famous 
relationship for the heat-transfer coefficient: 



Nju 



= 0.725 



n^^P(P-0^)h^-g 



1/4 



UD (T ^-T ) 
sat w 



( 2 . 1 ) 



In order to simulate and analyze a tube bundle, 
Eissenberg [Ref. 11] stated the following additional 
assumptions: 

8. Condensate drains as a laminar sheet from a tube on 
to the tube directly underneath in such a way that 
velocity and temperature gradients are not lost in 
the fall between tubes. 

9, The saturation temperature and the tube-wall 
temperature are constant for all tubes in the bank. 

Jakob [Ref. 19] extended the Nusselt analysis for 
filmwise condensation heat transfer on a vertical in-line 
row of horizontal tubes as shown in Figure la. 

The above-mentioned assumptions were combined with the 
assumption of constant temperature drop across the 
condensate film for all the tubes, and the average 
coefficient for a vertical row of N tubes was predicted to 
be: 
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UN D (T 
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Figure la. Figure lb. 

Idealized Condensation More Realistic Picture of 

on Banks of Tubes Condensation on Banks of Tubes 
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Upon, dividing equation (2.2) by equation (2.1), the Nusselt 
theory can be expressed as: 




(2.3) 



Equation (2.3) can be also expressed in terms of the 
local coefficient for the N-th tube: 



In reality, condensate does not drop off in a continuous 
laminar sheet, but drops off instead by discrete droplets of 
liquid, as shown in Figure lb, depending upon the surface 
tension of the condensate. These droplets create ripples in 
the condensate film, and thereby decrease the performance 
degradation due to inundation. Based on his research, Kern 
[Ref. 20] proposed a less conservative relationship: 




(2.4) 



\ - 1/6 
H7 = 



(2.5) 



or, in terms of the local coefficient for the N-th tube 



Chen [Ref. 21] considered the following conditions: 

1. the momentum gain of the falling condensate between 
tubes f and 

2, the condensation of vapor on the condensate between 
tubes , 

and concluded that: 



\u 



(1 + 0.2C(N-1) ) ( 



1-0.68? + 0.02?f. 
1+0. 95? - 0.15??^ 




(2.7) 



where; 



C = 



kAT 




? = 







and 




The above approximate expression, due to Chen, is valid for 
most ordinary applications. 

Experimental work doen by Eissenberg [Ref. 11], in order 
to investigate the effects of steam velocity, condensate 
inundation, and noncondensable gases on the heat-transfer 
coefficient, revealed that condensate does not always drain 
onto tubes aligned vertically, but can be diverted sideways, 
caused by local, vapor-flow conditions. The condensate thus 
follows a staggered path as shown in Figure 2. 

Eissenberg, making the assumption that the flow is 
dominated by gravity, stated in this side-drainage model, 
that condensate strikes the lower tubes on their sides 
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Figure 2. Droplet Path Through a Tube Bundle with Side Drainage 
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rather than their tops. Therefore, the inundation effects 
influence the condensate flow only on the lower half of the 
tubes, which transfer less heat than the upper half. Based 
on the above-mentioned model, Eissenberg obtained the 
following formula: 

\ -1/4 

^ = 0.60 + 0.42 N (2.8) 

Much experimental research has been conducted for 
studying the effect of condensate inundation. In general, 
the obtained data are highly scattered as shown in Figure 3. 
Berman [Ref, 22] made a comprehensive comparison of filmwise 
condensation data on bundles of horizontal tubes, and he 
concluded that the wide variation in expermental data for 
tube bundle inundation is caused by the following variables: 

1. bundle geometry (in-line or staggered), 

2. tube spacing, 

3. type of condensing fluid, 

4. operating pressure. 

5. heat flux, and 

6. local vapor velocity. 

In addition, noncondensable gases, and insufficient steam 
for lower tubes can cause such scattering of data. 
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Inundation Studies 



Nobbs [Ref, 18] used one active tube in a dummy tube 
bundle. He simulated additional condensate by using three 
porous tubes. Based on his results, he concluded the 
following; 

1, Vapor velocity increases the condensate heat- 
transfer coefficient on both inundated and 
uninundated tubes in a given tube bundle, 

2, The effect of inundation is to reduce the heat- 
transfer coefficient, 

3, The condensate drainage path is often not vertically 
downwards but in a diagonal direction. This can 
result in tubes receiving different amounts of 
inundation, 

Marto and Nunn [Ref, 23] made a comprehensive survey of 
the effects of noncondensable gases. They stated that these 
can be classified into one of two categories; 

1, the introduction of an additional local thermal 
resistance due to the propagation of noncondensable 
gases towards a condensing tube surface under the 
influence of a gradient in partial pressure, and 

2, the cumulative effect of gas blanketing where uneven 
rates of condensation in a condenser bundle eventu- 
ally lead to regions where tubes are inoperative in 
a condensing role. 

Experiments done recently by Cunningham [Ref, 15] revealed 
that noncondensable gases have a less effect on the perform- 
ance of finned tubes compared to smooth tubes. Based on the 
above discussion, it is clear that the performance degrada- 
tion due to noncondensables must be taken into consideration 
in realistic designs of tube bundles. 
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Another very important factor in the performance of a 
condenser is the effect of vapor velocity. As noted 
earlier, the vapor shear plays a beneficial role, Berman 
and Tumanov [Ref. 24] conducted experiments on a single 
horizontal tube placed in a bank of uncooled neighboring 
tubes. For vapor in vertical downflow, they found a 
relation between the vapor Reynolds number, and the heat 
flux as: 



= 1 + 9.5 X 10"^ (Rev) 
N 



11 • 8 / •\J~Nu 



with the restriction, that 



11.8/-Jnu 

(Rev) <50 



Eissenberg [Ref. 25] has stated that in designing exper- 
imental bundles, the combined effects of inundation and 
vapor shear are very important. In fact, the use of narrow 
condenser bundles to experimentally study vapor shear and 
inundation effects is preferred to simulate large conden- 
sers, However, small narrow tube bundles can create errors, 
due to the following causes: 

1. Wall flow: condensate drainage may reach side 

walls ; 

2. Dummy tubes: condensate inundation may disperse; 

3. Vapor lanes: vapor may bypass along the side walls; 

and 

4. Noncondensable gases: they can affect condensation 

even at low concentration in the bulk stream, 
particularly if steam is recycled. 
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III. EXPERIMENTAL FACILITY 



A. TEST FACILITY 

The test facility shown in Figures 4 and 5, was designed 
and built by Morrison [Ref. 13] and modified by Noftz 
[Ref. 12] to simulate an active tube column having up to 30 
tubes in incranents of five tubes deep (i.e., five, ten, 
fifteen, etc.). Some elements of the original test facility 
were modified to allow for more efficient operation of the 
facility. 

A detailed description of the components used in the 
test facility is given in Reference 12. Only a short 
description of these components will be found in this 
report. Particular attention, however, will be focused on 
the experimental tubes. Calibration procedures for 
components requiring calibration are outlined by Reilly 
[Ref. 7] . 

B. STEAM SYSTEM 

The steam system shown in Figure 6 was modified from 
Noftz 's initial design. From the house supply line, steam 
flows via a 19 mm O.D. stainless steel line through a steam 
supply valve (MS-3) to a cast iron steam separator. The 
steam continues through the system past two Nupro bellows 
valves, which were used in conjunction with the supply valve 
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Figure 4. Front View of Test Facility 
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Figure 5. Rear View of Test Facility 
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Figure 6. Schematic of Steam System 



to regulate the steam supply pressure. From these valves, 
the steam flows into the test condenser diffuser. 

The steam supply pressure was monitored by a pressure 
gage just downstream of the steam supply valve; also a 
compound gage just after the Nupro valves was used to 
monitor the pressure in this line. 

The operator had no control over the state point, 
quality or noncondensable gas content since house steam was 
used. However, especially at nights and during weekends, 
the state point of the steam at the inlet to the test 
condenser was found to be nearly constant. 

C. TEST CONDENSER 

The test condenser shown in Figure 7 was unchanged from 
Noftz's initial design. Steam enters via the top, passes 
through the transition piece, the vortex annihilator, and 
finally through the diffuser. The dimensions of the test 
condenser were 305 mm X 305 mm X 79 mm and it was made of 
stainless steel. These dimensions allowed for a maximum of 
twenty-seven 16 mm O.D. tubes arranged in an in-line config- 
uration of three columns of nine tubes each. For this 
experiment, the in-line configuration was used with a middle 
column of five active tubes flanked on either side by a 
column of five dummy tubes. Just on top of the upper active 
tube, a perforated, distilled water supply tube was posi- 
tioned, and flanked on both sides by dummy tubes. In order 
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NOTE- ALL COMPONENTS DRAWN TO SCALE 



Figure 7. Sketch of Test Condenser 
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to conduct the experiments, a square, in-line arrangement of 
tubes, with a pitch-to-diameter ratio of 1,5 was used. 

A vertical slot (Figure 7) along the centerline of each 
condenser end plate was used for active and perforated tube 
installation. The tubes were positioned using nylon tube 
sheets that were attached to the exterior of the condenser 
end plates. 

To minimize heat losses, and also to prevent leaks from 
the tube sheets, each test condenser side was provided with 
a nylon tube sheet one-inch thick with six holes 
(S/D = 1,5), having about 0.5 mm tube clearance, which 
allowed the tubes to be easily slid into the test condenser 
through the tube sheets. The exterior side of each tube 
sheet had grooves to support 0- rings. The aluminum tube 
sheets had six holes, having a tube clearance of about 
0,5 mm and were used as sealing plates. The diffuser, 
exhaust plenum, transition piece, vortex annihilator, and 
the exhaust piping which are shown in Figures 8 and 9, were 
insulated with rubber insulation. 

A viewing window allowed viewing of the condensation 
process, A double-walled glass window was used, and heated 
air was fed through the clearance between the two glasses to 
eliminate fogging on the inner glass. 



36 




Figure 8. Details of Transition Piece and Vortex Annihilator 
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EXHAUST (STEAM AND CONDENSATE) 

f 




TO TEST CONDENSER HOTWELL 



Figure 9. Details of Exhaust and Condensate Piping from 
Exhaust Plenum 
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D. TEST CONDENSER TUBES 



Four kinds of active tubes were tested in this experi- 
ment. The tubes were manufactured by the Wolverine Division 
of Universal Oil Products, and all were made of titanium. 

The first kind of tubes, shown in Figure 10 were smooth 
titanium tubes of 16 mm O.D. with a 1,65 mm wall thickness. 
The second kind of tubes, shown in Figure 11, were single- 
start, helically-corrugated tubes, designated as Low 
Pressure Drop (LPD) , with 16 mm O.D. and a 1.65 mm wall 
thickness. The third kind of tubes, shown in Figure 12, 
were also single-start, helical ly-corrigated LPD tubes, 
wrapped with titanium wire of 1.58 mm O.D, The fourth kind 
of tubes, shown in Figure 13, were smooth tubes, wrapped 
with 1.58 mm O.D. titanium wire on the same pitch as the 
roped Wolverine tubes. In order to wrap the titanium wire, 
the tubes were attached into a lathe, and the wire was 
welded at the start of the helical groove. Then, under a 
constant tension of 4 kg, the wire was wrapped manually by 
turning the chuck of the lathe. Finally, the free end of 
the titanium wire was welded at the end of the helical 
groove. 

The 660-mm-long active tubes were connected to separate 
cooling water supply and discharge lines. The dummy tubes 
flanking the active tubes were made of 16 mm O.D. stainless 
steel. These tubes served to direct the steam flow so as to 
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Figure 10. Photograph of Smooth Titanium Tube 
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Figure 11. Photograph of Roped Titanium Tube 
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Figure 12. Photograph of Wire Wrapped Roped Titanium Tube 
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Figure 13. Photograph of Wire Wrapped Smooth Titanium Tube 






simulate actual conditions in a condenser. Cooling water 
was not supplied to these tubes and they did not penetrate 
the test condenser end plates. Special characteristics of 
the Wolverine tubes are listed in Table I, 

E. PERFORATED TUBE 

The perforated tube water supply system is shown in 
Figure 14, This system consisted of a perforated tube 
(located above the uppermost active tube ) r a water heater 
which served as a supply tank, a rotameter in order to 
control the amoung of water supplied to the perforated tube, 
a pump driven by a 1/2-HP electric motor, a condensate pump 
and associated piping system and valves. The active length 
of the copper-nickel, perforated tube was 305 mm, which was 
identical to the length of the test condenser. Supply water 
entered one end of this tube; the other end was sealed off. 

The supply tank was unchanged from Noftz's initial 
design, except for the addition of a thermocouple in order 
to have a direct indication of the exact temperature of the 
water supplied to the perforated tube. The water, heated to 
the temperature of the condensate leaving the bottom tube, 
was fed to the perforated tube by a 1/2-HP, electric-motor- 
driven pump, via a rotameter and a recirculation valve. The 
flow rate to the perforated tube was controlled by using the 
rotometer and the valve provided in the heater-water 
recirculation line. Also, another modification was made to 
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TABLE I 



DESIGN CHARACTERISTICS OF WOLVERINE 
CORODENSE TUBES (TYPE LPD) 



Outside diameter (in) 


: 


5/8 


Wall thickness (in) 


• 


0.035 


Outside diameter (ft) 


: 


0.0520 


Heat Transfer Surface and Ratio 






a. Outside A* (ft^/ft) 




0.163 


* * 

b. Outside to inside A /A. 

o' 1 




1.141 


Inside diameter (ft) 




0.0435 


2 

Cross section for flow A (ft ) 

cs 




0.00163 


Number of groove starts 




one +0 


Pitch (in) 




0.300 


Groove radius (in) 




1/32 nominal 


Depth-transition length 
(next to plain section) (in) 




1-3/4 max 
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Figure 14. Schematic of Perforated Tube Water Supply System 





feed the condensate collected in the hotwell, via the 
condensate pump, to the supply tank (heater) in order to 
facilitate the inundation of up to 30 tubes under atmos- 
pheric conditions. 

The main feature of the perforated tube was to inundate 
with condensate from above, in order to simulate a tube 
column of more than five active tubes. When the wrapped 
Wolverine tubes and the smooth wrapped tubes were tested, 
the perforated tube was also wrapped with titanium wire, by 
using the above-mentioned technique. 

F. CONDENSATE SYSTEM 

The condensate system shown in Figure 15 was modified 
from Noftz's initial design, and consisted of the test 
condenser and hotwell, the condensate pump, piping and 
valves . 

The test condenser hotwell collected the condensate 
produced by the test tubes. With valve C-1 closed, the 
condensate mass flow rate from the test condenser could be 
measured using a stop watch. Opening the valve, the 
condensate was fed, via the condensate pump, to the supply 
tank for the perforated tube, or dumped into the building's 
drainage system. The condensate line connecting the test 
condenser and the test condenser hotwell were insulated 
using rubber insulation. 
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Figure 15. Schematic of Condensate System 




■ k- 




G. COOLING-WATER SYSTEM 



The cooling-water system was a partially-closed system 
as shown in Figure 14, and it was unchanged from Noftz's 
initial design. House water was used for the test facility. 
Cooling water was stored in a 1,2 meter, cubical, plexiglass 
supply tank, and was pumped by a 5-HP, electrically-driven 
pump via 51-mm-O.D, plastic piping to a manifold. Five 
rotameters were attached to the manifold to measure the flow 
rate through each active tube. Five regulating valves were 
used to obtain any water velocity between 0 and 5 m/s within 
the tubes. The cooling water passing through the rotameters 
was fed via 16 mm O.D. stainless steel tubing, and tygon 
tubing to the active tubes. The total tube lengths were 
long enough (over 2 meters) to ensure hydraulically 
fully-developed flow, and no swirling into the tubes. After 
flowing through the tubes, the cooling water was passed 
through mixing chambers, where the temperature profile was 
destroyed, to facilitate the steady measurement of the 
outlet cooling water temperatures. After the mixing 
chambers, the cooling water was collected in the supply 
tank, A 7, 5-HP, electrically-driven pump was used to pump 
water from the supply tank, through a filter, to a cooling 
tower in order to minimize temperature rise at the water 
inlet. 
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The cooling tower was located outside the building and 
was composed of four truck radiators across which air was 
blown by means of a fan. The entire system, consisting of 
the heat exchanger and the fan was enclosed in a wooden 
structure with louvered openings to provide enough 
ventilation. 

The tygon tubing was secured to the inlet and outlet 
sections of the active tubes by means of hose clamps. The 
outlet sections, including the mixing chambers, were 
insulated with rubber insulation. 

H. INSTRUMENTATION 
1 , Flow Rates 

a. Foulton rotameters were used to measure the flow 
rate of cooling water for each active tube. Starting from 
the top active tube to the bottom one, the rotameters were 
calibrated giving 100% maximum flow rates of 66.9+1, 

72,6+1, 71.5+1, 72.8+1, and 73.3+1 kg/min. 

b. The perforated tube water rotameter was 
calibrated giving a 100 percent volumetric flow rate of 
969±30 ml/min. 

c. All rotameters were calibrated using the 
procedures noted in Appendix A of Reference 7, 
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2. Temperature 



Stainless-steel sheathed, copper-constantan 
thermocouples were used as the primary temperature 
monitoring devices. Three thermocouples were used to 
measure the inlet cooling water temperature; and twelve 
thermocouples were utilized to measure the outlet cooling 
water temperature. 

Additionally, two thermocouples were used to measure 
the steam saturation temperature; one thermocouple was used 
to monitor the condensate temperature, and one was used to 
measure the vapor temperature. Table II lists the locations 
monitored, 

A Gulton Industries, West 20, 0-500° F temper- 
ature controller was used to regulate the temperature of the 
perforated tube supply water. The controller had a 
manufacturer-stated accuracy of +0.5% of the span (about 
±1.25° F) . It was impossible to obtain the precise 
temperature for the perforated tube water, because of the 
low response of the heating system. 

During inundation runs, the condensate temperature 
rose with the increasing number of tubes. For example, the 
condensate temperature after the 5th tube was 94.6° C 
and it increased to 97.65° C after the 30th tube. 
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3 . Pressure 



Several different types of pressure measurement 
devices were used in this facility. They were: a Bourdon 

tube pressure gage which was used to measure the steam 
supply pressure, located downstream of the steam supply 
valve; an absolute pressure transducer which was used to 
measure the test condenser pressure; and, a compound 
pressure gage was also used to measure the pressure drop 
downstream of the Nupro valves. The pressure transducer was 
calibrated against a mercury manometer, 

4 . Data Collection and Display 

A Hewlett-Packard HP-3054A Automatic Data 
Acquisition System, with a HP-9826 computer and a HP-2671G 
printer were used to record and display the thermocouple and 
pressure transducer readings. The temperatures were 
recorded in degrees Celsius, while the pressure was recorded 
in volts which was then converted, using the calibration 
curve, into mm Hg absolute. The pressure transducer was 
assigned to channel 019 of the HP Automatic Data Acquisition 
System, while the thermocouples were assigned channels as 
indicated in Table II. 
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TABLE II 



CHANNEL NUMBERS FOR COPPER-CONSTANTAN THERMOCOUPLES 



Location 



Tci »1 

T^i #3 

■"oi *5 

CO 

T^^ #1 

CO 

T^_ #1 
CO 

T ^ #2 

CO 

T^^ #2 

CO 

T„ #3 

CO 

T ^ #3 

CO 

T^^ #4 

CO 

CO 

T ^ #4 

CO 

T ^ #5 

CO 

CO 



sat 



sat 



cond 



'vap 



Channel 

000 

001 

002 

003 

004 

005 

006 

007 

008 

009 

010 

011 

012 

013 

014 

015 

016 

017 

018 
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IV. PROCEDURES 



A. INSTALLATION AND OPERATING PROCEDURES 

1 * Pi:eparation_of_CQnd ensei^T u bes 

Prior to installation, the titanium tubes were 
cleaned using a chemical cleaning method [Ref. 26]. The 
steps in this cleaning process are as follows; 

a. Swab the tube surface with acetone to remove grease. 

b. Using a test tube brush, brush the inside surface of 
the tube with a 50% sulfuric acid solution in order 
to remove any oxides- Also, apply this solution 

to the outside surface of the tube. 

c. Rinse the inside and outside of the tube with tap 
water . 

d. Apply, using a brush, a 50% solution of sodium 
hydroxide mixed with an equal amount of ethyl 
alcohol, at the boiling temperature (about 

85° C) to the outside surface of the tube, 

e. Rinse the tube with tap water. 

f. Rinse thoroughly with distilled water. 

Prior to any run- the condenser tubes had to be prepared to 
ensure filmwise condensation- Exterior and interior 
surfaces were cleaned to ensure proper wetting character- 
istics. Also, the tubes were cleaned by running steam at 
atmospheric pressure through the test condenser for about 
twenty minutes without the cooling water running through the 
tubes . 
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It was also found that, when drop-wise condensation 
occurred, rinsing the tubes using the perforated-tube water 
was sufficient to restore film-wise condensation. 

2 . System Operation and Steady-State Conditions 

Complete operating instructions are listed in 
Appendix A, The steady-state condition was reached about 
three hours after initial system light-off, and about 
fifteen minutes after changes to the cooling water or 
perforated tube supply water flow. 

When a steady-state condition was reached, the runs 
were made. The duration of each run was approximately one 
minute. For each condition, for example, for tubes 11 
through 15, five consecutive runs were made and average 
values were computed. 

The following data were taken automatically by the 
data-acquisition system; 

a. the thermocouple readings, and 

b. the pressure transducer reading. 

Also, the following data were read into the computer, 
through the keyboard, for each run: 

a. the setting of each rotameter, and 

b. the test condenser hotwell levels. 
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B. DATA REDtlCTION PROCEDURES 

1 . Overall Heat-Transfer Coeffi c ient 

The heat-transfer rate to the cooling water is given 

by: 

Q = mC (T -T . ) f 4 1) 

^ PW CO Cl 

The heat-transfer rate can also be found from the 
overall heat-transfer coefficient by: 



Q = U A • LMTD (4.2) 



where 



LMTD 




-T . ) 
Cl 



sat 



-T 



sat 



-T 




(T^-T^^) 
S CO 



(4.3) 



After combining equations (4.1), (4.2) and (4.3), it 
is found that 



U 



o 




T ^-T . \ 
sat Cl \ 

T ^-T 
sat CO / 



(4.4) 



2. Inside Heat -T ransfer Coefficient 

The heat-transfer coefficient on the inside is 
calculated from the Sieder-Tate relationship as described in 
Holman [Ref . 27 ] : 
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r 



1/3 ^ji_j 0.14 



y ' 
w 



(4.5) 



In the above equation. is referred to as the 



Sieder-Tate coefficient. The remainder of the riqht-hand 
side of the above equation. 



is referred to as the Sieder-Tate parameter, X. 

For thermally and hydrodynamically developed flow in 
tubeSf the Sieder-Tate coefficient equals 0.027. However, 
when short tubes are considered, as in the case of this 
experiment, fully-developed conditions are not attained. 
Therefore, the value of must be found experimentally. 



the Sieder-Tate coefficient. The Wilson plot was developed 
in 1915 [Ref. 281. It is merely a plot of l/U^ versus the 
inverse of the Sieder-Tate parameter (which is proportional 
to the inverse of cooling water velocity raised to the 0.8 
power) . The reasoning behind the Wilson plot is shown in 
the following development. 

Consider the equation for the overall heat-transf er 
coefficient: 



(R 



e 



0.8 p 1/3 y 0.14 
• r ( — ) 



The Wilson plot was used to arrive at the value of 
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( 4 . 6 ) 




1 



o 



or 




h 



1 



o 



( 4 . 7 ) 



where 




D 

o 



R 



w 



For the Wilson-plot method to be successful, h^ 
must be kept constant. When steam condenses on the shell 
side, the above condition can be achieved only if the heat 
flux (q") is kept constant. 

As required for the Wilson plot, when the cooling 
water velocity increases- h^ increases. increases, 
and finally q" increases - 
Consider 



decrease LMTD by lowering the steam saturation temperature. 
This requires tr ial-and-error setting of condenser pressure 
which is a very difficult task. 



q” = UqLMTD 



( 4 . 8 ) 



In order to keep q" constant, one must therefore 
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To avoid this difficulty and yet arrive at a Sieder 
Tate constant not affected by the varying q", a modified 
Wilson-plot method was used. This method was developed by 
Wanniarachchi [Ref. 291, and the required steps are listed 
below: 

1. Assume = 0.027 

2. Calculate 



LMTD 



«,n 




-T 



ci 



sat 



-T 



sat 



-T 




3. Calculate 

Q = me (T - T^ . ) 

P CO Cl 

4 . Calculate 



A ♦ LMTD 



5A. For the first data point only, do the following 

a. Assume 



Cf = O = 1 

U 



b. Calculate the Sieder-Tate parameter 



X = R p (iL) 



c. Calculate 



h. 

1 




R 

e 



0.8 p 1/3 
r 



(^) 
y ' 
w 



0.14 
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d. Calculate 




e. Calculate 



f. Repeat steps a through e until C^ 
assumed in step a approximately equals 

calculated in step e. 

g. Calculate 

1-1 1 

h U h. A. w 

o oil 

h. Assign = Q 



NOTE: The second subscript of h on the left-hand side of 



equation (4.9) refers to the first data point. 

5B. For all data points e xcept for first data 
point, do the following: 

a. Calculate 



h. 

1 



o 



_1 

U 



R 



w h 



b,N 



1 




where 




b. Calculate 
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c 



Calculate the Sieder-Tate parameter 




r 



(_L) 0-14 



6. Repeat steps 2 through 5 for all data points. 



7. Plot ^ versus 



1 

X 



NOTE: It is more reasonable to plot h. versus X. 

However, 1/h. versus 1/X was piotted in order to 
be consistent!: with the original Wil son-plot method 



versus X. 



( ^ versus ^ ) • 



o 



8. Obtain the slope (m) . of the least-squares-f it . 
straight line. 

9. Calculate the Sieder-Tate constant. 



10. Repeat steps 2 through 9 until the assumed and 
calculated values are approximately equal. 

A listing of the computer program used in this 



method can be found in Appendix D. 

3 . Outside Heat-Transfer Coefficient 

The heat-transfer coefficient on the outside is 
calculated using the following steps: 

1. Calculate the average bulk water temperature. 



D. 



C. 

1 



1 
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2. Evaluate the thermophysical properties based on 
the average bulk water temperature. 

3. Calculate the cooling water velocity. 



V = 
w 



m 

pA, 



4. Calculate the water-side Reynolds number. 



R = 
e 



P V D. 
w w 1 



w 



5. Calculate the heat transferred to the cooling 
water . 



Q = m (T - T . ) C 

CO cx pw 



6. Calculate the heat flux, 
Q 



q” = 



TT* D • L 
o 



7. Calculate the Nusselt coefficient using the 
formula: 



h„ = 0.651 
Nu 




(4.10) 



NOTE: The above formula was employed since no direct 

measurement of tube wall temperature was made. 

To compare the condensate film temperature, as 

required for equation (4.9), an iterative scheme was used as 

outlined below: 
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NOTE: 



a. Assume = T 

f sat 

b. Evaluate the relevant thermophysical properties, 
which are included in equation (4.10) . 

c. Calculate the Nusselt coefficient using 
equation (4.10) - 

d. Evaluate T_ using the formula: 

X # c 



T 



ffC 



sat 




5 



e. Repeat steps a through e until assumed 
in step a approximately equals calcu- 

lated in step e. 

Calculate the inside heat-transfer coefficient using 
the formula: 




(4.11) 



In order to determine the C^, the average inner 
wall temperature must be known. As noted earlier, 
this tCTiperature is not known directly and it must be 
found iteratively, as described below: 

a. Assume a correction factor C^ (say 1) . 

b. Calculate the h^ using equation (4.11) . 

c. Calculate the tube wall temperature using the 
formula : 




d. 



Calculate a new correction factor, at the 
evaluated tube wall temperature 



C 



f 




0.14 
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e. Repeat steps a through d until assumed 
in step a approximately equals calcu- 
lated in step d. 

9. Calculate the log-mean- temperature difference, 
(LMTD) . 



LMTD 




(4.12) 



10. Calculate the overall heat-transfer coefficient 
using the formula: 

^ LMTD (4.13) 



11. Calculate the outside heat-transfer coefficient 
using the equation: 



h 




D 



Dihi 



R 



w 



(4.14) 



NOTE: The validity of equation (4.13) is based on the 

assumption of negligible water-side fouling 
resistance and the resistance due to noncondensables. 

12. Calculate the normalized, local, outside heat- 
transfer coefficient, 

13. Calculate the normalized average, local outside 
heat-transfer coefficient (h^/h^) . 
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C. DATA- REDUCT ION PROGRAM 



A computer program was utilized to analyze the rav; data. 
The program was in BASIC language and was run on an HP-9826 
computer system, A peripheral plotter was used to plot the 
results. The program is presented in Appendix D. 
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V. RESULTS AND DISCUSSION 



A. SIEDER-TATE COEFFICIENTS FOR SMOOTH AND ROPED TUBES 

The modified Wilson Plots for smooth and roped tubes are 
shown in Fiqures 16 and 17. Of particular interest is the 
very small intercept, which represents the estimated fouling 
factor- For example, the estimated fouling factor for data 
run STSD-11 (for smooth tubes) is 1.37X10“® m^^K/W 
based on the inside area. This value is only three percent 
of the typical value (4.4X10”® m^^K/W) used for 
design purposes. This very small value supports the initial 
assumption of negligible fouling factor for the success of 
the modified Wilson-Plot method. 

Table III is a summary of the Sieder-Tate coefficients 
calculated for smooth and roped tubes. The average Sieder- 
Tate coefficient was calculated to be 0,029±0.001 for 
smooth tubes and 0.061+0.002 for roped tubes. Thus, the 
roped tubes have a Sieder-Tate coefficient 2.1 times greater 
than that for the smooth tubes. This increase is mainly due 
to the increased surface area, turbulence and swirl effects. 

The Sieder-Tate coefficient derived for the smooth tubes 
is 7.4 percent greater than the value (0.027) published in 
the original, generalized correlation. This increase can be 
easily explained by the short condensing tubes used in this 
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Figure 16. Modified Wilson-Plot for Smooth Tubes (Run STSD-11) 
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Figure 17. Modified Wilson-Plot for Roped Tubes (Run RTSD-11) 



TABLE III 



SUMMARY OF SIEDER-TATE COEFFICIENTS 
FOR SMOOTH AND ROPED TUBES 



File Name 


Tube Type 


Sieder-Tate 

Coefficient 


STSD-2 


Smooth 


0.0287 


STSD-10 


Smooth 


0.0291 


STSD-11 


Smooth 


0.0296 


RTSD-3 


Roped 


0.0591 


RTSD-4 


Roped 


0.0589 


RTSD-6 


Roped 


0.0621 


RTSD-7 


Roped 


0.0627 
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experiment. Even though the flow condition was hydro- 
dynamically fully-developed, it was thermally developing 
throughout the condensing length, resulting in a greater 
Sieder-Tate coefficient. 

It is worth noting that the Sieder-Tate coefficient 
derived by using the original, Wilson-Plot method (i.e., 
plotting 1/U^ versus 1/X) , consistently gave values about 
10% greater than the values obtained using the modified, 
Wilson-Plot method. 

B. SMOOTH TUBES 

Figure 18 shows the variation of the normalized, local, 
condensing coefficient for five tubes. The data points lie 
about 40 percent above the curve predicted by the Nusselt 
theory. This close agreement was considered to be an 
indication of the proper opepation of the test apparatus and 
the data reduction procedures. 

As discussed earlier in Chapter II, the Nusselt theory 
for a tube bundle is based on a number of basic assumptions. 
For example, in reality, the condensate does not fall as a 
continuous sheet as assumed for the Nusselt theory, but it 
forms drops before leaving the tube. This phenomenon has a 
tendency to result in a larger condensing coefficient than 
that predicted by the Nusselt theory. 
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Figure 18. Variation of Local Condensing Coefficient with Tube Number 
(Run STNWNI-1) 



Figure 19 shows the normalized, average condensing 
coefficient for five tubes. These data points show a much 
smoother trend than that for local values. Further, these 
data points are closer to the Eissenberg correlation than 
the Nusselt prediction. 

Figure 20 shows the variation of the normalized, local, 
condensing coefficient for 30 tubes. The data points lie up 
to 35 percent above the curve predicted by the Nusselt 
theory- The outside heat-transfer coefficielt for the first 
tube was ten percent smaller than the value predicted by the 
Nusselt theory. 

Figure 21 displays the normalized, average condensing 
coefficient for 30 tubes. Again, these data points show a 
smoother trend than that for local values. These points lie 
between the Eissenberg and Nusselt predictions. 

Figure 22 shows the least-squares-f it curve for the data 
points for smooth tubes. The exponent derived is -0.154, 
and it is in close agreement with the value of -0.14 derived 
by Noftz [Ref. 12] . 

C. ROPED TUBES 

Figure 23 displays the variation of the normalized, 
local, condensing coefficient for five tubes. The data 
points lie up to 15 percent above the curve representing the 
Nusselt theory- 
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Figure 19 . Variation of Average Condensing Coefficient with Tube Number 
(Run STNWNI-1) 
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Figure 20. Variation of Local Condensing Coefficient with Tube Number 
(Run STNWI-1) 
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Figure 21. Variation of Average Condensing Coefficient with Tube Number 
(Run STNWI-1) 
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Figure 22. Least-Squares-Curve Fit for Data {Run STNW-1 ) 
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Figure 23. Variation of Local Condensing Coefficient with Tube Number 
(Run RTNWNI-3) 



Figure 24 shows the normalized, average, condensing 
coefficient for five tubes. These points show a trend very 
close to the Eissenberg relationship. 

Figure 25 represents the plot of the normalized, local, 
condensing coefficient under inundation conditions up to 30 
tubes . 

Figure 26 shows the normalized, average, condensing 
coefficient for 30 tubes. These data points lie about 
midway between the Eissenberg and Nusselt predictions. 

Figure 27 shows the least-squares-f it curve for the data 
points for roped tubes. The exponent calculated is -0.183. 
Since the interest is in the large tube bundles, the curve 
fit was generated only for tubes 11 through 30. This is 
evident from the figure as the curve fit shows poor agree- 
ment when compared to the data points for the first ten 
tubes . 

D. ROPED TUBES WRAPPED WITH WIRE 

The outside, heat-transfer coefficient for the first 
tube was 11.3 percent greater than the value predicted by 
the Nusselt theory. This increase is in agreement with the 
manufacturer's claim that the special geometry of the roped 
tubes has a thinning effect on the condensate film, result- 
ing in a larger condensing coefficient. 
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Figure 24. Variation of Average Condensing Coefficient with Tube Number 
(Run RTNWNI-3) 
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Figure 25. Variation of Local Condensing Coefficient with Tube Number 
(Run RTNWNI-3) 
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Figure 26. Variation of Average Condensing Coefficient with Tube Number 
(Run RTNWI-1) 
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Figure 27. Least-Squares-Curve Fit for Data (Run RTNWI-1 ) 



Figure 28 shows the variation of the normalized, local, 
condensing coefficient for five tubes. The data points lie 
up to 63 percent above the curve predicted by the Nusselt 
theory. 

Figure 29 displays the normalized, average, condensing 
coefficient. The data points are well above the Eissenberg 
correlation. 

Figure 30 shows the normalized, local, condensing 
coefficient for a bundle of 30 tubes. The data points are 
scattered within the limits of uncertainty, as predicted by 
the error analysis (see Appendix C) . 

Figure 31 shows the variation of the data points 
representing the normalized, average, condensing coefficient 
under inundation up to 30 tubes. The data points lie up to 
100 percent above the curve predicted by the Nusselt theory. 

Figure 32 represents the least-squares-f it curve, which 
has a derived exponent of -0.039. The curve fit is in good 
agreement with all the data points. 

E. SMOOTH TUBES WRAPPED WITH WIRE 

The outside heat-transfer coefficient for the first tube 
was up to six percent greater than the value predicted by 
the Nusselt theory. This increase is caused by the thinning 
effect on the condensate film, resulting from the surface- 
tension forces acting toward the wrapped wire. 
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Figure 28. Variation of Local Condensing Coefficient with Tube Number 
(Run RTWNI-3) 
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Figure 29. Variation of Average Condensing Coefficient with Tube Number 
(Run RTWNI-3) 
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Figure 30. Variation of Local Condensing Coefficient with Tube Number 
(Run RTWI-3) 
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Figure 31. Variation of Average Condensing Coefficient with Tube Number 
(Run RTWI-1) 
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Figure 32. Least-Squares-Curve Fit for Data (Run RTWI-1 ) 




Figure 33 shows the variation of the normalized, local, 
condensing coefficient for five tubes. The data points lie 
up to 76 percent above the curve predicted by the Nusselt 
theory- 

Figure 34 shows the normalized, average condensing 
coefficient for five tubes. These data points are well 
above the Eissenberg correlation. 

Figure 35 shows the variation of the normalized, local, 
condensing coefficient for 30 tubes. The data points lie up 
to 107 percent above the curve predicted by the Nusselt 
theory. 

Figure 36 displays the normalized, average, condensing 
coefficient for 30 tubes. The data points lie above the 
curve representing the Eissenberg correlation. 

Figure 37 shows that the curve fit is in good agreement 
with all the data points. The derived exponent is -0.037. 

F. OBSERVATIONS 

1 . Smooth Tubes 

During all runs, complete film-wise condensation was 
observed without any visible evidence of drop-wise conden- 
sation. This was mainly achieved by the tube cleaning 
procedures used. 

A slow rate of condensate droplet migration, from 
cooling water outlet end to the inlet end, was observed at 
the bottom of each active tube. This problem was minimized 
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Figure 33. Variation of Local Condensing Coefficient with Tube Number 
(Run STWNI-2) 
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Figure 34. Variation of Average Condensing Coefficient with Tube Number 
(Run STWNI-2) 
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Figure 35. Variation of Local Condensing Coefficient with Tube Number 
(Run STWI-1) 
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Figure 36. Variation of Average Condensing Coefficient with Tube Number 
(Run STWI-1) 
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Figure 37. Least-Squares-Curve Fit for Data (Run STWI-1) 



TABLE IV 



SUMMARY OF RUNS WITH INUNDATION 



File Name 


Tube Type 


External Wire 
Wrapped or Not 


STNWI-1 


Smooth 


No 


RTNWI-1 


Roped 


No 


RTNWI-2 


Roped 


No 


RTNWI-3 


Roped 


No 


RTNWI-4 


Roped 


No 


RTNWI-5 


Roped 


No 


RTWI-1 


Roped 


Yes 


RTWI-2 


Roped 


Yes 


STWI-1 


Smooth 


Yes 


STWI-2 


Smooth 


Yes 
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TABLE V 



SUMMARY OF RUNS WITHOUT INUNDATION 



File Name 



Tube Type 



External Wire 
Wrapped or Not 



STNWEV-2 


Smooth 


No 


STNWNI-1 


Smooth 


No 


STNWNI-2 


Smooth 


No 


STNWNI-3 


Smooth 


No 


RTNWNI-3 


Roped 


No 


RTWNI-1 


Roped 


Yes 


RTWNI-2 


Roped 


Yes 


RTWNI-3 


Roped 


Yes 


STWNI-2 


Smooth 


Yes 


STWNI-4 


Smooth 


Yes 


STWNI-5 


Smooth 


Yes 
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RESULTS FOR RUNS WITHOUT INUNDATION 
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TABLE VI (continued) 



VD 


CN 


VO 


CM 


in 




m 




• 


• 


• 


• 


rH 


CM 


rH 


CM 



rH 


ro 


in 


VD 


in 


iH 




O 


o 




in 


in 


a\ 


00 


a\ 


a\ 


• 


• 


• 


• 


o 


o 


o 


o 





VO 




CM 


o 


<J> 


rH 


CM 






rH 


ro 


00 


00 


a\ 


a\ 


• 


• 


• 


• 


o 


o 


o 


o 



w 


cn 


cn 


cn 


<D 


<D 


0) 


(U 




>H 




>H 




1 


T 


in 

1 


1 

H 




1 

H 








CO 




cn 



98 



RESULTS FOR RUNS WITH INUNDATION UP TO 30 TUBES 
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COMPARISON OF HEAT-TRANSFER COEFFICIENTS FOR TUBE #1 IN THE BUNDLE 
(WITHOUT INUNDATION) AT 1.56 m/s COOLING WATER VELOCITY 
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TABLE VIII (continued) 
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COMPARISON OF HEAT-TRANSFER COEFFICIENT FOR TUBE WITH INUNDATION 
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TABLE IX (continued) 
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TABLE X 



COMPARISON OF FOR UNINUNDATED TUBE RUNS 



Tube Type 


External Wire 
Wrap or not 




Smooth 


No 


0.7986 


Roped 


No 


0.9019 


Roped 


Yes 


0.9651 


Smooth 


Yes 


0.9339 
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TABLE XI 



COMPARISON OF h FOR INUNDATION TUBE RUNS 

Xu 



Tube Type 


External Wire 
Wrap or not 




Smooth 


No 


0.59 


Roped 


No 


0.63 


Roped 


Yes 


0.86 


Smooth 


Yes 


0.86 
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TABLE XII 



EXPONENTS OF THE LEAST-SQUARES-FIT 



File Name 


Tube Type 


External Wire 
Wrap or not 


Exponent 


STNWI-l 


Smooth 


No 


0.154 


RTNWI-1 


Roped 


No 


0.183 


-2 


Roped 


No 


0.191 


-3 


Roped 


No 


0.179 


-4 


Roped 


No 


0.185 


RTWI-1 


Roped 


Yes 


0.039 


-2 


Roped 


Yes 


0.039 


STWI-1 


Smooth 


Yes 


0.037 


-2 


Smooth 


Yes 


0.056 
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after leveling the tubes, by adjusting the leveling nuts on 
the test condenser support bracket- When the condensation 
or the inundation rate increased, it was observed that the 
drops were formed at more sites along the tubes, but the 
droplet size for each tube was nearly the same. 

2 . Roped Tubes 

Again, during the experiments, there was no evidence 

of drop-wise condensation. It was observed that the falling 

drops were coalesced on the lower surface of the tube, 

especially in the space between the two successive grooves. 

and then rivulets were formed and fell on the tube below. 

Under condensate inundation, the phenomenon of the rivulet 

formation was more intense. It is also worth noting that 

the droplet-formation frequency was higher near the cooling 

water inlet end. This can be easily explained by the 

larger, local, Sieder-Tate coefficient and by the larger 

temperature difference, T - , at the inlet compar- 

ed to T - T , at the outlet end. 
sat CO 

The condensate was formed in the space between the 
successive grooves wrapped with wire. Then due to the 
surface tension forces, the condensate was drawn to the base 
of the wire, and drops were formed at the bottom surface of 
the tubes. 
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Under inundation conditions, the above-mentioned 



drop formation and movement were more intense. It is worth 
noting that no splashing was observed during the experiment. 

4. Smooth. Tubes Wrapped with Wire 

It was observed that the condensate forms between 
two successive, helically-wrapped wires, and was then drawn 
immediately towards the space between the wire and the tube 
surface. Rivulets were drawn from the base of the wire to 
the next tube below. It was observed that there was no drop 
migration along the tube. 

Under condensate inundation, the above mentioned 
drop formation was more intensive, but again no splashing 
was observed. 
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VI. 



1. The average, outside, heat-transfer coefficient for 
30 smooth tubes was 59 percent of the Nusselt coefficient 
calculated for the first tube in the bank. 

2. The average, outside, heat-transfer coefficient for 
30 smooth tubes wrapped with wire was 86 percent of the 
Nusselt coefficient calculated for the first tube in the 
bank , 

3. The average, outside, heat-transfer coefficient for 
30 roped tubes was 63 percent of the Nusselt coefficient 
calculated for the first tube in the bank. 

4. The average, outside, heat-transfer coefficient for 
30 roped tubes, wrapped with wire was 86 percent of the 
Nusselt coefficient calculated for the first tube in the 
bank , 

5. The Sieder-Tate coefficient for roped tubes was 2.1 
times greater than that for smooth tubes. 

6. Wire wrapping considerably improves the average 
condensing coefficient for both smooth and roped tubes in 
tube bundles. 

7. Of all cases investigated in this study, roped tubes 
with wire wrap would be the best candidate for designing 
compact condensers. However, the water side pumping power, 
which was not investigated in this work, could be consider- 
ably higher than that for the smooth- tube case. 
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VII. RECOMMENDATIONS 



A. TEST APPARATUS MODIFICATIONS 

The following test apparatus modifications are 
considered advisable: 

1. Redesign the test condenser flanges to allow a large 
sealing area, and also to allow for the possibility 
of further experiments with different S/D ratios. 

2. Redesign the existing test condenser hotwell to 
allow more reliable and convenient measurement of 
condensate flow rate. 

3- Install a larger heater for the perforated-tube 
water supply tank to facilitate the simulation of 
larger tube bundles. Also, modify the temperature 
control system for the pe rf or ated-tube water supply 
tank, to allow for more rapid heating and cooling to 
reduce long delays between runs. 

4. Install a glass window on the rear side of the test 
condenser, in order for the operator to have a view 
on both sides of the test tubes. 

5. Make the system vacuum tight so that data can be 
taken at vacuum conditions'. 



B. ADDITIONAL TESTS 

The following additional tests would be important in 
this continued investigation. 

1. Conduct tests with enhanced tubes manufactured by 
Yorkshire Imperial Metals, Ltd. 

2. Take movies of the condensation process so that 
further conclusions could be drawn with regard to 
condensate drop phenomena and their relationship to 
condenser performance. 

3. Conduct tests varying the pitch of the wrapped wire 
on smooth titanium tubes and determine the optimum 
wire pitch. 
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4 . 



Withdraw gas samples from the test condenser and 
analyze on the gas chromatograph to determine the 
effect of noncondensable gases on the condenser 
performance . 

5, Conduct tests to investigate the effect of vapor 
velocity on the heat-transfer coefficient. 
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APPENDIX A 



OPERATING PROCEDURES 



A. INITIAL PROCEDURES 

1. Energize the main circuit breaker located in Power 
Panel P-2 on the wall to the right of the test 
apparatus. 

2. Energize the circuit breaker on the left side of the 
old control board by pressing the ON button. 



Energize the following switches 
panel : 


in the control 


a. 


#1 - Perforated tube water 
pump) . 


supply tank (feed 


b. 


#2 - Outlets. 




c. 


#3 - Hot water heater. 




d. 


#4 - Condensate pump. 




e. 


#6 - Cooling tower. 




f . 


#7 - Cooling water pump. 





4. Ensure all test apparatus valves are closed. 

5. Fill the perforated tube condensate supply tank with 
distilled water. Set the temperature controller for 
the perforated tube water supply tank at about 

95° Q, fully open the recirculation valve. P-1, 
and start the pump to begin heating the water. The 
controller will have to be reset to the proper 
supply temperature once steady-state conditions are 
obtained . 

6. Fill the cooling water supply tank. This can be 
done by backfilling with valves CK-1 and CW-4 open. 

7. Energize the data acquisition system. 
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B. OPERATION 

1 . House Steam 

a. Open the main supply valve. 

b. Open valve MS-3 until the pressure gage 
indicates the desired steam supply pressure. 

c. Fully open MSD-1 to drain any condensate. 

d. Open valves MS-4 and MS-5 until the desired 
steam supply pressure is obtained and re-adjust 
MS-3 as necessary. 

2. Condensate System, To collect the condensate in the 
test condenser hotwell, operate the system with 
valve C-1 closed. After a test run is completedr 
open valve C-1 to drain the condensate by opening 
valve C-2 to the bilges, or operating the condensate 
pump in order to fill the perforated tube supply 
tank , 

3 . Cooling-Water System 

a- Open valyes CW-1 . CW-2 and CW-3 . 

b. Ensure valves CKl-1 and CW-4 are closed. 

c. Energize the two cooling water pumps. 

d. Open valves CW-5 . CW-6 . CW-7 . CW-8 and CW-9 to 
obtain the desired cooling rates, 

4. Perforated Tube Water Supply System 

a. Once steady-state conditions are achieved, reset 
the temperature controller to the proper 
inundation temperature. 

b. Adjust the rotameter to the required flow rate 
for each run. The supply tank recirculation 
valve may have to be adjusted to achieve the 
desired flow rate, but should never be fully 
closed to avoid damaging the pump and to ensure 
uniform water temperature in the tank. 

c. Refill the supply tank as required, by using the 
existing piping filling system, by operating the 
condensate pump; or filling the supply tank with 
distilled water filling line. 
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5. Miscellaneous. 



To maintain a clear test condenser window, open 
valve A-2 and then energize and adjust the air 
heater power supply. When securing, always turn off 
the power supply first and allow the air heater to 
cool before securing valve A-2 . 



C. SECURING THE TEST APPARATUS 

1. Secure the steam valves MS-5, MS-4, MS-3 and the 
main supply valve- 

2. Secure the air compressor. 

3. Secure the pe rf or ated-tube water supply system by 
securing the pump, temperature controller, and 
valves P-1 and P-4. Drain the supply tank by open- 
ing valve P-3 (if desired). 

4. Secure the test condenser viewing window air heater 
as prescribed above. 

5. Secure the data acquisition system. 

6. Allow the test condenser to cool down for about 30 
minutes- then secure the cooling water pumps and 
close valves CS-1, CW-2 , CW-3 , CW-5 . CW-6 , CW-7 , 
CW-8, CW-9 and C-4. 

7. Drain the test condenser hotwell. 

8. Secure all circuit breakers. 

9. Drain the cooling water system piping and rotameters 
by opening valve CW-10 and leaving open the cooling 
water rotameter supply valves - 

10. Drain the cooling water supply tank by opening the 
drain valve via the remote operating rod. 

11. Ensure all valves are secured. 
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APPENDIX B 



SAMPLE CALCULATIONS 



A. RUN STNWNI-1 



A sample calculation is performed in this section to 
illustrate the solution procedure used in the data reduction 
program [Ref 29] . 



The STNWNI-1 run was selected to perform this analysis; 
INPUT PARAMETERS 



File 

Pressure condition 
Inundation condition 
Month, date and time 
Run number 

Tube number 1 2 

Inlet temp (°C) 28.53 28.54 

Outlet temp (OC) 32.55 32.20 

Saturation temperature 
Degree of superheat 
Condensate temperature 
Static pressure 



STNWNI-1 
Atmospheric 
5 tube s 

05:11:10:08:50 

1 

3 4 5 

28.55 28.54 28.53 

32.17 32.05 31.83 

100.24 (°C) 

1.88 (OC) 

92.41 (OC) 

767.33 (mm Hg) 



The following calculations are further limited only to 



the first tube. 



1. Determination of Average Bulk Temperature 



T,(l) = T^,-(l) + T^^(l) xO.5 

b cl CO 

Tj^(l) = (28.53 + 32.55) 0.5 

T, (1) = 30.540c 
o 
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2 . Thermophysical Properties 



P = 5.329 
r 



P = 995.2 kg/m' 



u = 798 X 10"^ N.s/m^ 



C = 4.178 Kj/kg.K 

pw J' 3 



k = 619 X 10 
w 



-3 



W/m-K 



NOTE: All properties are calculated at the average bulk water 

temperature from Table A. 6 p. 782 [Ref 30] 

3 . Cooling Water Mass Flow Rate 
m = 14.22 Kg/min = 0.237 Kg/s 

4 . Determination of Cooling Water Velocity 



V = — = 
w PA. 



V = 
w 



(14.22) 60 



.-4, 



V 



w 



(995.2) (1.56 X 10 ) 

= 1.53 m/s 



Determination of Reynolds Number 



^e = 



p V D. 
w w 1 

P 

w 



(995.2) (1.53) (0.141) 



798 X 10 



-6 
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= 26,904 



6, Determination of Heat Transfer 



Q = m • (T - T . ) •€ 

CO Cl pw 



Q = (14.22) (^) (32.55 - 28.53) (4,178) 



Q = 3,980.5 W 
Determination of Heat Flux 



q" = 



7T.D .L 
o 



„„ ^ 3,980.5 

^ ^-(0.015875) (0.305) 



q" = 261,682.2 

8. Determination of Nusselt Coefficient 






= 0.651 



kf^. Pf^.hfg. (9.81) 



lif.Do.q" 



1/3 



Assume T, = T . 

f sat 



= 100.24 C 



= 957.8 kg/m' 



h. = 2.2564 X lO"' J/kg 

fg 



= 682 X lO”^ W/m-K 
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NOTE; 



9. 



10 . 



NOTE ; 



-6 N 

u = 281 X 10 ^ 

2 

m 

All properties are calculated at the Saturation 
Temperature. From Table A. 6 p. 782 [Ref 30] 






= 0.651 



■ 32 1 1/3 

(682x10"^) • (957.8) (2256.4x10^) (9.81) 



(281xl0“^) (0.015975) (261682.2 



= 11492.6 W/m^K 



Determination of T 



f ,c 



= ^sat - hj,^ 



T. = 100.24 - 0.5 

f,c 11492.6 

T. =88. 850C 
f ,c 



Thermophysical Properties 



k = 673 X 10"^ W/mK 



p = 963 Kg/m^ 

U = 326 X 10"® N-s/m^ 
h^g = 2.289.5 X 10^ J/kg 

All properties are calculated at the film temperature. 
From Table A. 6 p. 782 [Ref 30] 
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11. Determination of Nusselt Coefficient 






= 0.651 



k 3 P 2h 
^f ’^f ^fg-9.81 
Uf.Do.q" 



1/3 



^Nu =0-651 



3 2 11/3 

(673 X 10 * (963) (2289.5 x 10 ) (9.81) 



(326 X 10“®) (0.015875) (261685.4) 



= 10,885.3 W/m .K 



12. Determination of T 



f,c 



T 



f /C 



’sat 



q” 

Nu 



0.5 



T 



f ,c 



100.24 



261685.4 

10885.3 



T. = 88.220C 
f / c 



13 . Determination of Logarithmic Meal Temperature 
Difference (LMTD) 



LMTD 



T 



CO 



T . 
Cl 



In 



( 



^sat~^ci \ 
"^sat '^co/ 



LMTD 



32.55 - 28.53 
,100. 24 - 28. 53^ 
^100.24 - 32.55^ 



LMTD = 69.68°C 



14. Determination of Overall Heat - Transfer Coefficient 



U = —3 — 

O LMTD 

261685.4 
% 69.68 
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u = 

o 


w 

^ 3755.5 

m^K 



15. Determination of Inside Heat-Transfer Coefficient 



Assume C£ = 1.1 





= 0.029 


h. = 
1 


!e.C r p 0-333 

D^ ^i*^e ’^r *^f 


h. = 
1 


0*0141 ^0-^29) (26904°*®) (5.829°*®®®®) -1.1 


II 

•H 


9303.2 W/m®K 



16 . Determination of Inner Wall Temperature 



T = 
w 


T + ^ 

b h. D. 


T = 
w 


,, , (261685.4) (0.015875) 

JU.04 + (9303.2) (0.0141) 


T = 
w 


62.20C 



17. Determination of at the Average Wall Temperature 



y (62.20C) = 453 x lO”® N-s/m^ 
w 

18. Determination of Correction Factor 



Cfc 


f-w 1 0.14 

L'‘„ (65-59-)J 


°^fc 


r 798 X 10'®1 
L453 X 10"°J 


^fc 


= 1.08 
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19. Determination of Outside Heat-Transfer Coefficient 



1 



U 



o 




R 



w 



h 



1 



1 



0.015875 



0.000042925 



3755.5 



(0.1014) (9303.2) 



h = 9772.3 W/m^K 



B. RUN STSD-11 

A sample calculation is performed in this section to 
illustrate the solution procedure used in the modified Wilson 
plot program [Ref 29] . 

The DP-11 run was selected to perform this analysis 
INPUT PARAMETERS 

File DP-11 

Month, date and time 04:01:13:57:20 

Data point #1 

Steam Saturation Temperature ' 100.01°C 

Inlet Temperature 24.42°C 

Outlet Temperature 30.75°C 

Flowmeter Reading 10% 

1 . Determiantion of Average Bulk Water Temperature 




T, = (24.42 + 30.75) 0.5 = 27.59°C 
b 



2 . Thermophysical Properties 



P 



r 



5.83 



m 
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c. 



pw 



4.180 Kj/Kg.K 



P = 857 X 10"® N-s/m^ 



3. 



4 . 



5 . 



k = 0.613 W/m-K 
w 

All properties are calculated at the average bulk water 
temperature from Table A. 6 p. 782 [Ref. 30] 

Cooling Water Mass Flow Rate 

= 6.36 kg/min 

Determination of Cooling Water Velocity 



^ ^ (6.36) 60 

^ (997) (1.56 X lO""^ 

V =0.68 m/s 
w 

Determination of Reynolds Number 



R 



P V D. 
w w 1 



e 



y 

w 



R 



e 



(997) (0.68) (0.0141) 

857 X 10“^ 



R^ = 11,154 



Determination of Heat Transfer 




Q = (6.36) (;^) (30.75 - 24.42) (4180) 



Q = 2,805 W 
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Determination of Logarithmic Mean Temperature 
Difference (LMTD) 



LMTD 



A n 



/ "^sat '^ci 
V'^sat '^co 



) 



LMTD 



30.75 - 24.42 
/lOO.Ol - 24.42\ 
UOO.Ol - 30.75/ 



LMTD = 72.380C 



8. Determination of Overall Heat - Transfer Coefficient 



U = i 

tt.d -L-LMTD 
o 



U = 



2805 



o 7 t(0. 015875) (0.305) (72.38) 

= 2,547.715 -Jf— 
m K 

9 . Determination of Sieder - Tate Parameter 



^ -0.8.P -0.3333 

e r 



X = (11,154"°-®) (5.83"°'^^^^) 

X = 0.0003212 

10 . Determination of Inside Heat - Transfer Coefficient, 
based on assumed Sieder - Tate Coefficient 



Assiame C^ = 1.1 and 
C^ = 0.03 



h. 

X 



^c.,^o.a.,^o.333c^ 
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613 X 10 
0.141 



) ( 1 . 1 ) 



11. 



12 . 



13. 



14. 



-3 



h. = 
1 



(0.03) (11,154®’^) (5.83°*^^^ 



W 



h. = 4,463.4 _ 

^ m^’K 

Determination of Average Wall Temperature 

T - T = Q 

W b 1T.D; .L*h. 

1 1 



T, . T, = 2805 

w 7r(0.0141) (0.305) (4,453.4) 

T, - T. = 46.510C 
w b 



T = 46.51 + T, 
w b 



T = 46.51 + 27.59 = 74°C 
w 



Obtain at the Average Inner Wall Temperature 



U (T ) = 375 X lO”^ N-s/m^ 
w 



Determination of 




C 



fc 



0.14 



^ ^ / 857 X 10 

\375 X 10 

C. = 1.12267 
fc 

Iterate for 
Assume Cj = 1.12 and 
= 0.03 
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I 




h = ^ C R “-®P “-”33.0 
"i D. ‘'e 
1 



h. 

1 



613 X 10 



-3 



0.0141 



(0.03) (11,154°-^) (5.83°-^^^^) (1.12) 



h. 

1 



= 4,546.97 W/m K 



15 . Determination of Wall Thermal Resistance, Based on 
the Outside Diameter 



R = 
w 



(D -D. ) D 
O 1 o 

k . (D +D. ) 
m o 1 



„ _ (0.015875 - 0.0141) (0.015875) 

~ (21.9) (0.015875 + 0.0141) 



R = 0.000042925 
w W 



Assume R^ = 0 
^ / 



16. Determination of Outside Heat - Transfer Coefficient 



h = 



U 



R - 

w D . h . 
1 1 



h = j 

2,547.715 

h = 10273.8 - 
m 

Set Q = Q = 
o 



1 

- 0.000042925 

W 

2.K 

2805 W 



0.015875 

(0.0141) (4,463.4) 



17. Determination of Actual Sieder-Tate Parameter 



X = 



X 




_ 0.0003212 

1.12267 
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X = 0.0002861 



INPUT PARAMETERS 



File STSD-11 

Data point #2 

Steam saturation temperature 100.03°C 

Inlet temperature 24.390C 

Outlet temperature 29.460C 

Rotameter setting 15% 



18 . Determination of Average Bulk Water Temperature 



T 



b 



(T . + T ) X 0.5 

Cl CO 



T, = (24.39 + 29.46) 0.5 
b 



Tj^ = 26.920C 



19 . Thermophysical Properties 

P = 5,823 
r 



u = 996 Kg/m^ 

Cp^ = 4.180 Kj/KgK 

u = 855 X 10“^ N.S/m^ 

k = 613 X lO"^ W/mK 
w 



20 . Cooling Water Mass Flow Rate 
M^ = 9.72 Kg/min 

22 . Determination of Cooling Water Velocity 
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22 . 



23. 



24. 



25. 



= 



(9.72) 60 



-4 



(996) (1.56 X 10 ’) 



V = 1.04 m/s 
w ' 



Determination of R, 



R = 
e 



P V D. 
w w 1 



^ (996) (1.04) (0.0141) 



855 X 10 



-6 



R^ = 17,082 
e 



Determination of Heat - Transfer 



Q = M- (T - T . ) C 

f CO Cl pw 



Q = (9.72) (|!^) (29.46 - 23.39) (4180) 

Q = 3,433.2 W 

Determination of Logarithmic Mean Difference 



Temperature (LMTP) 



LMTD = 



T - T . 

CO Cl 



iln 



/ "^sat *^ci \ 
V^sat-^co / 



LMTD = 



29.46 - 24.39 
. /100.03 - 24.39\ 
\100.03 - 29.46/ 



LMTD = 73.08°C 

Determination of Overall Heat-Transfer Coefficient 



tt.d .L-LMTD 
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• (0.015875) (oTsOS) (73.08) 
= 3,088.4 W/m^K 



26. Determination of Sieder-Tate Parameter 



X 




p -0.3333 
r 



X = (17,082"°-®) (5.823"°*^^^^) 

X = 0.000228 

27. Determination of ^ 







) 




f_i_ 

L3088.4 



1 

10273 



^ 2805 ^ 

\3,433.2/ 



1/3 

-0.00004295 



0.015875 

0.0141 



i = 0.0002095 JLS 



28 . Determination of Average Wall Temperature 



T = T + = 

w b r -D. .L*h, 



T - 26 92 + 3433.2 

iw - + Ti. (0-015875) (0.305) (4,773.26) 



T = 74.20OC 
w 



29 



Determination of at the Average Wall Temperature 



p 



w 






375.2 X 10”^ N-s/m^ 
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30. Determination of 



C 



fc 



_!iw_ 



855 X 10 



-6 



0.14 



'fc 



375.2 X 10 



-6 



= 1.1222 



31. Determination of Actual Sieder-Tate Parameter 



X 



X 



C 



fc 



^ _ 0.000228 
^ " 1.12459 

X = 0.0002027 
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APPENDIX C 



UNCERTAINTY ANALYSIS 



The general form of Kline and McClintock [Ref. 31] "second 
order" equation is used to compute the probable uncertainty in 
the results. For some resultant, R, which is a function of 
primary variables ■ t the probable uncertainty in 

R, 6 R is given by: 



6 R = 



<ilr ‘1: “ 2 ' 



^ , 0R \2 

• ^ex- ^^n> 



1/2 



(C.l) 



^1 - -n 

where <SX^ , 6 X 2 ,..., is the possible uncertainty in each of 

the measured variables. 



A. UNCERTAINTY IN THE COOLING WATER VELOCITY 



V = 



m 



w P A^ 



Applying equation (C.l) the following equation results: 



6 V 



w 



V 



w 



4,2 , ^ ,^,2] 



1/2 



(C.2) 



The following uncertainties were assigned to the variables: 

6 m = + 0 . 01 kg/s 
3 

6 p = +3 kg/m 



5A^ = +0.0001 m 
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B. UNCERTAINTY IN THE REYNOLD'S NUMBER 









w 1 



w 



The probable uncertainty is given by: 



R 



o 6V ^ 6D. ^ X n 

(|^)^ + (-^)^ + (-5^)^ + 

p n 



1/2 



(C.3) 



The following uncertainties were assigned to the variables 
6p = +3 kg/m^ 



6D^ = +0.0001 m 

6|a = +8x10 ^ N*s/m^ 

6V^ = from equation (C.2) 



C. UNCERTAINTY IN HEAT TRANSFER 



Q = m(T -T . )C 
^ CO Cl pw 



The probable uncertainty is given by: 



1/2 



iQ = 
Q 



• r >p 

, 6m , 2 , , CO 



O 6T . - 6C - 

,)2 , ci _,2 ^ ,_^)2 



T -T . 

CO Cl 



(C.4 ) 



CO Cl CO Cl pw 

The following uncertainties were assigned to the variables: 
6m = +0.01 kg/s 



6T = +0.025 °C 

CO — 

ST . = +0.025 °C 

Cl — 

SC = +8 J/Kg*C 
pw — ^ 
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D. UNCERTAINTY OF THE HEAT-FLUX 



q" = 



ttD L 
o 



The probable uncertainty is given by; 



1/3 



6q'' ^ 1 



^ q ' ' d 



0)2 ^ ^ 6^)2 



(C.5) 



The following uncertainties were assigned to the variables 



6D = +0.0001 m 
o — 

6L = +0.0001 m 

6Q = as found from equation (C.4) 



E. UNCERTAINTY OF h„ 

Nu 

1/3 



The probable uncertainty is given by: 



h„ = 0.651 
Nu 



■, 3 2. 

MfDoq 



h 

h 



Nu 

Nu 



0.651 



, 4 f ,2 . , 2«£,2 . , 4^,2 H. ',^ 2,2 , 
kf 3 p 3h^g 3g 



f .2 



3l^^ 



+ 



+ ( 



6D 

1d^ 



o , 2 



. 1/2 



) + (■ 



3q 



tt) 



(C.6) 



The following uncertainties were assigned to the variables: 
6k^ = +0.0012 w/m*k 



2 

6p = ±3 kg/m 
6h^g = +0.48 J/kg 
6Hf = +8x10"^ N*S/m^ 
= +0.0001 m 
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6g = +0.001 m/s^ 

6q" = as found from equation (C.5) 



F. UNCERTAINTY OF T. . , 

film. 



T. = T ^ 
fc sat 



Nu 



0.5 



The probable uncertainty is given by: 




V^sat '2 is:. 2 ,!V)2' 

T..^'^ ^ q"^ ^ h. ■* 



1/2 



sat 



Nu 



(C.7) 



The following uncertainties were assigned to the variables: 

5T ^ = +0.025 °C 
sat — 

6q"= as found from equation (C.5) 

6h,, = as found from equation (C.6) 

Nu 



G. UNCERTAINTY OF OVERALL HEAT-TRANSFER COEFFICIENT 
a" 

U = 

o LMTD 

The probable uncertainty is given by: 



1/2 



U. 



I 5 ( LMTD ) , 2 
^ q"^ ^ LMTD ' 



(C.8) 



The following uncertainties were assigned to the variables: 
6q = as found from equation (C.5) 

6 (LMTD) = as found from equation (C.9) 



133 



H 



UNCERTAINTY FOR LOGARITHMIC MEAN TEMPERATURE DIFFERENCE 
( LMTD ) 



LMTD = 



T -T . 

CO Cl 

1 -/ s Cl , 

S CO 



The probable uncertainty is given by: 



5 LMTD 
LMTD 



6T (T .-T ) 

s Cl CO 



T -T . T^-T^^ T -T . 

^ S Cl S CO , , s Cl,, 

ln(= — = — ) 



T -T 
s CO 



6T 



Cl 



(T -T . ) T -T . 

S Cl , , S Cl,/ 

s; — ) 



T -T 

S CO 



6t 



^ 1/2 



CO 



(T -T ) T -T . 

. S CO , , S Cl ,/ 

In ' rn ' 



T -T 

S CO 



(C.9) 



The following uncertainties were assigned to the variables: 
6Tg = +0 .025 °C 
6T . = +0.025 °C 

Cl — 

6T = +0.025 °C 

CO — 



I. UNCERTAINTY OF INSIDE HEAT-TRANSFER COEFFICIENT 



h. 



1 



kw . ^ rjO-Sn 0.333, u,0.14 
D^ ^i^e ^r Vw 
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The probable uncertainty is given by: 



6h. 



6k „ 6D. 

+ (-^) + ( 



w 



D. 

1 



0.86R _ 0.3336P 

^ 



R 



0.146(n/n ) ~ 

+ ( rrrr: 



1/2 






m 



- <SC. - 
-)^ + (-^)^ + 
^i 



(C. 10) 



The following uncertainties were applied to the variables: 

6k = +0.0012 w/m*k 
w — 

6D. = +0.0001 m 
1 — 

6C^ = +0.0001 

6Rg = as found from equation (C.3) 

6Pj. = +0.17 

6(|a/|j ) = 8x10 ^ Ns/m^ 
w 



J. UNCERTAINTY IN TEMPERATURE DIFFERENCE 



AT = 



h. D. 
1 1 



6AT 

AT 



(^,2 + (!!) i )2 ( 1^.2 

^ q" ^ ^ h. ^ D 

^ 1 o 



<Sd. - 



D. 

1 



1/2 



(C.ll) 



The following uncertainties were assigned to the variables 

6q" = as found from equation (C.5) 

6h^ = as found from equation (C.IO) 

6D = +0.0001 m 
o — 



6D^ = +0.0001 m 
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K. UNCERTAINTY IN OUTSIDE HEAT-TRANSFER COEFFICIENT 



h = T- 
o 1^ 

U. 



D. 



- R 



Dihi 



w 



h 



+ 



.U' 



6U 



1 D 

(1 - R - ) 

' U D . h . 

O 11 

6R. ^ 



w 



D 

— R — ^ 

D . h . 

o 11 



6h. 



Dihi 



) (■ 



D^ 

D . h . 

o 11 



2 ^ 



1/2 



(C.12) 



The following uncertainties were assigned to the variables: 
6U^ = as found from equation (C.8) 

6hi = as found from equation (C.IO) 



6R = +0.00001 m -k/w 
w — 



L. UNCERTAINTIES FOR THE NORMALIZED LOCAL HEAT-TRANSFER 

COEFFICIENT h.,/h . 

N' 1 

This ratio is simply the heat transfer coefficient of a 
given tube, N, divided by that of the first tube, i.e., for 
the fifth tube, N=5 and: 

% 




136 












4 



An application of equation (C.l) results in the following 
equation ; 



6(h„/hi> 




(C.13) 



Note; Equation (C.13) is valid only for N7/2. For example; 



M. UNCERTAINTY IN THE NORMALIZED AVERAGE HEAT-TRANSFER 
COEFFICIENT, 

The normalized average heat-transfer coefficient is obtained 
for the Nth tube by taking the average of the heat-transfer 
coefficients of the first N tubes and dividing this by the 
heat-transfer coefficient of the first tube; 



Applying equation (C.l) to the above, the following equation 
results : 



/h ) 





hjq (h^ + h2 



+ . . . + hj^ ) N 




2 



(yhi) 

(V^l) 



(C.14) 



i=l 



where N = the tube number . For example ; 
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For Run STNWNI-1 



= 1.54 + 0.11 "'/s 

Re = 27546 + 2121 

Q = 3934 + 169 W 

q" = 267830 + 11516 / 2 

^ ^ m 



tr, = 10802 + 164 "/ 2, 
— 'm.Jc 



LMTD = 69.7 + 0.6 °C 

U = 3843.7 + 165 2, 

0 — '^m*k 

w 

h. = 9381 + 122 V 2, 

1 — m.k 




11361.4 + 
= 0.8891 



W 

1249 / 2 

+ 0.14 




0.9446 + 0.12 



APPENDIX D 



COMPUTER PROGRAMS 



A. DATA REDUCTION PROGRAM 



1000 ! 
1010 ! 
1020 
1030 
1040 
1050! 
lOGO! 
1 070! 
1080! 
1090 
1100 
1110 
1120 ! 
1 130! 
1140! 
1 150 
1160 
1170 ! 
1180! 
1190! 
1200 
1210 
1220 ! 
1230! 
1240 
1250 
1260 
1270 
1280 
1290 
1300! 
1310! 
1320 
1330! 
1340! 
1350 
1360 
1370! 



FILE MAME: DRP 
REVISED: May 20. 1983 

COM /Cl/ C(7) 

DIM Tci<2) .Tco(4.2) .Ti(4) .Mf t(4) .Vu-<4> .Ho(4) 

DIM To(4) ,Ts(l ) ,Tb(4) .R3(4) ,R4(4) .$3(4) .S4<4) 

ASSIGN COEFFICIENTS FOR THE 3-TH ORDER 
POLYNOMIAL FOR TYPE-T (COPPER-CONSTANTAN) 
THERMOCOUPLES 

DATA 0. 10086091 .25727 .94369 . -767345 . 3295 . 78025595 . 81 
DATA -9247486589. 6. 97688E+1 1 , -2 . 661 92E+1 3 .3 . 94078E+1 4 
READ C<*) 

assign Full-scale flgu rates th.rouch the 5 
FLOW METERS (kQ/mn) 

DATA 66.36.73.35.72.44.72.52,72. 24 
READ Mft(<^) 

ASSIGN SIEDER-TATE COEFFICIENT AND EXPONENT 

FOR REYNOLDS NUMBER 

Ci=.029 

Ex=.8 



ASSIGN GEOMETRIC VARIABLES 



Di=.0141 
Do=.0I5875 
Ktn=21 .9 
L=.305 
Nc=3.3333 
Pt=1 .5 



Inner diameter <m) 

Outer diameter <m) 

Thermal conductivity of titanium (W/m-K) 
Condensing length <m) 

Number of unit cells across condenser width 
Transverse tube pi tch- co-diameter ratio 



COMPUTE THE MINIMUM STEAM FLOW AREA IN THE TEST CONDENSER (m'2) 
Amf=Nc*Do*'(Pt-PI/(4»Pt ) )*L 



COMPUTE INSIDE AREA AND WALL RESISTANCE 
Ai-PI*Di"2/4 

Rid=Do»LQG (Do/Di)/(2»Ktm) 



1380 PRINTER IS 701 
1390 CLEAR 709 
1400 BEEP 

1410 INPUT "ENTER MONTH. DATE. AND TIME < MM : DD ; HH :MM ; SS)" . T imeS 
1420 OUTPUT 709:"TD":TimeS 
1430 BEE? 

1440 INPUT "ENTER THE INPUT MODE < 1 =3054A ,2=FILE)'Mm 
1450 IF Im=2 THEN 
1460 BEEP 

1470 INPUT "ENTER THE NAME OF THE EXISTING BATA FILE" .OlddataS 
1480 PRINT USING "10X.""This analysis was performed for data stored 
10A":01ddataS 

1490 ASSIGN ?File2 TO OlddataS 
1500 END IF 
1510 IF Im=l THEN 
1520 BEEP 

1530 INPUT "GIVE A NAME FOR THE DATA FILE TO BE CREATED" . NewdataS 
1540 CREATE BOAT NewdataS. 20 
1550 ASSIGN jFilel TO NewaataS 



in 



file " 
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1560 END IF 
1570 BEEP 

1580 INPUT "GIVE A NAME FOR THE OUTPUT FILE".File_oatS 
1590 BEEP 

1600 INPUT "ENTER THE PRESSURE CONDITION U=ATM.2»VACUUM)" .Mp 
1610 IF Md= 1 THEN PRINT " Pressure condition: ATMOSPHERIC" 

1620 IF Mp=0 THEN PRINT " Pressure condition: VACUUM" 

1630 BEEP 

1640 INPUT "ENTER THE INUNDATION CONDITION (1=5 TUBES. 2=30 TUBES)". Mi 
1650 IF Mi=2 THEN PRINT " Inundation condition: 30 TUBES" 

1660 IF Mi=1 THEN PRINT " Inundation condition: 5 TUBES" 

1670 CREATE BOAT Fiie_outS.6 

1680 ASSIGN ?File3 TO File_outS 

1690 Ja=0 

1700 Nrun=0 

1710 FOR 1=0 TO 4 

1720 S3<I)=0. 

1730 S4ili=u. 

1740 NEXT I 
1750 Repeat: ! 

1760 Nrun=Nrun+I 
1770 OUTPUT 709:"TD" 

1780 ENTER 709;Time$ 

1790 PRINT " " 

1800 PRINT USING " 1 0X .""Month, date, and tine: "’M5A":TimeS 
1810 IF In=2 THEN Rdf 
1820 BEEP 

1830 INPUT "ENTER FLOW METER READINGS (AS PERCENTAGES)" .Fnl .Fn2 
1840 IF Nrun MOD 5=1 AND Mi =2 AND Hrun>5 THEN 
1850 BEEP 

I860 INPUT "ENTER FLOW RATE FOR POROUS TUBE (AS A PERCENT) ". Fpt 

1870 OUTPUT 'SFiieT:Fpt 

1880 Mpt=-8. 361613+10. 076742*Fpt 

1890 END IF 

1900 DISP "START COLLECTING CONDENSATE" 

1910 BEEP 
1920 WAIT 20 

1930 OUTPUT 709:"AR AFO AL19" 

1940 OUTPUT 722;"F1 R1 T1 Z1 FL1" 

1950! 

I960! READ INLET WATER TEMPERATURES 
1970! 

1980 FOR 1=0 TO 2 
1990 OUTPUT 709:"AS SA" 

2000 ENTER 722;Tci(I) 

2010 CALL Tvsv/(Tci(I>) 

2020 Tci ( I)=FNTemp( Tci ( I) . I) 

2030 NEXT I 
2040! 

2050! READ OUTLET WATER TEMPERATURES 
2060! 

2070 Ii=2 

2080 FOR 1=0 TO 4 

2090 IF 1=0 OR 1=3 THEN 

2100 Iu=2 

2110 ELSE 

2120 Iu=1 

2130 END IF 

2140 FOR J=0 TO lu 

2150 Ii-Ii+1 

2160 OUTPUT 709: "AS SA" 
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21 70 

2180 

2190 

2200 

2210 

2220 ! 

2230! 

2240! 

2250 

2260 

2270 

2280 

2290 

2300 

2310! 

2320! 

2330! 

2340 

2350 

2360 

2370 

2380! 

2390! 

2400! 

2410 

2420 

2430 

2440 

2450! 

2460! 

2470! 

2480 

2490 

2500! 

2510! 

2520! 

2530 

2540 

2550 

2560 

2570 

2580! 

2590! 

2600! 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720! 

2730! 

2740! 

2750 

2760 



ENTER 722:Tco(I.J) 

CALL Tvsv < Tco< I . J 1 ) 

Tco( I , J)=FNTemp( Tco<I,J) ,Ii) 

NEXT J 
NEXT I 

READ STEAM TEMPERATURES 

FOR 1=15 TO 16 
OUTPUT 709: "AS SA" 

ENTER 722:Ts(I-15) 

CALL Tvsv < Ts ( I” 1 5) > 

Ts( I - 15) =FNTenip( Ts< I~15) ,1) 

NEXT I 

READ CONDENSATE TEMPERATURE 

OUTPUT 709: "AS SA" 

ENTER 722:Tcon 
CALL Tvsv<Tcon) 

Tcon=FNTepp< Tcon. 1 7) 

READ VAPOR TEMPERATURE 

OUTPUT 709: "AS SA" 

ENTER 722 :Tv 
CALL Tvsv(Tv) 

Tv=FNTenp< Tv , 1 8) 

READ VAPOR PRESSURE 

OUTPUT 709: "AS SA" 

ENTER 722:P_volts 

COMPUTE AVERAGE WATER TEMPERATURES AT INLET 



Ti<0)=Tci<0) 

Ti(1 > = <Tci<0)<-Tci(1 ))*.5 
T i<2) = Tci< 1 > 

Ti(3>=<Tcl<1)+Tci(2))».5 

Ti(4)=Tci<2) 

COMPUTE AVERAGE WATER TEMPERATURES AT OUTLET 

FOR 1=0 TO 4 
IF 1=0 OR 1=3 THEN 

To(I)=(Tco(I .0)+Tco<I , l)+Tco(I,2) )*. 3333 
ELSE 

To(I)=<Tco(I.0)+Tco(I.1)>*.5 
END IF 
NEXT I 

Tsa=(Ts(0)+Ts(1>>*.5 
Pvap=FNPvsv(P_vol ts) 

Tsat = FNTvsp< Pvap) 

Dsup=Tv-Tsat 

READ INFORMATION FOR CONDENSATE FLOW RATE 
BEEP 

INPUT "ENTER INITIAL AND FINAL LEVELS IN HOT WELL 1".H1.H2 
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2770 Dh=H2-H1 

2780 IF Hran MOD 5=1 THEN Msun=0 
2790 Mf1=540.A836»Dh 

2800 Md1=.1f !-FMRhou<Tsat-IO>»l .OE-6/60 

2810 1 

2820 IF Mi=2 AND NmnOSO AND Nmn MOD 5-0 THEN 
2830 Mave = Ms(im/5 

2840 Se t = ( Mave-FNRhouC Tsat - 1 0 > / 1 0 *6+ . 03238) / . 042 1 32 
2850 END IF 
2860! 

2870 Rdf: ! 

2880! 

2890 PRINT USING "10X.""Ran number - •*" ,DD" ; Nrun 

2900 PRINT Tube ^ : 1 2 3 4 5" 

2910 IF Im=2 THEN 

2920 IF Nrun _M0D 5=1 AND Mi=2 AND Nrun>5 THEN ENTER s>File2;Fot 
293U lNTlk 'jr 1 ie2 : T 1 ( * ) . To( ~ > , Tsa. Toon , T v .Pvap , Teat , Dsup , Fn’ . Fiti 2 
2940 ENTER 5File2:HI.H2 
2950 END IF 

2960 IF Im=I THEN OUTPUT ?Fi le i : Ti ( » ) . To( * ) . Tsa. Tcon . Tv .Pvap . Tsat , Dsup .Fml ,Fm2 
2970 PRINT USING ’*1 OX . ""Inlet temp (Deg C) : .5( DDD . DD . 2X)" : T i < *) 

2980 PRINT USING " 1 OX . ""Outlet temp (Deg C) : "" ,5( DDD .DD . 2X) " : To( -) 

2990 PRINT USING "1 OX , ""Saturat ion temoerature = "".3D.DD."" (Deg C) ;Tsat 

3000 PRINT USING "1 OX . ""Degree of superheat = "".3D.DD."" (Deg C>”"";Dsup 

3010 PRINT USING "10X. ""Condensate temperature = "",3D.DD,“" (Deg C) :Tcon 

3020 PRINT USING "1 OX , ""Stat ic pressure = "".3D.DD."" (mm Hg) iPvap 

3030 IF Im=1 THEN OUTPUT ftFile1;H1.H2 
3040! 

3050! CALCULATE AVERAGE BULK TEMPERATURES 
3060! 

3070 FOR 1=0 TO 4 

3080 Tb(I)=(Ti(I)+To(I>>».5 

3090 NEXT I 

3100! 

3110 IF Hi=1 OR (Hi=2 AND Nrun<6) THEN As!=0. 

3120 IF Mi=2 AND Nrun>5 THEN S1=As1 

3130 S1=Asl 

3140 FOR J=0 TO 4 

3150 IF J=0 THEN Cuf=FmI 

3160 IF .J-1 THEN Cuf=Fra2 

3170 Mf=Mft(J)=Cwf/( 100=60) 

3180 Tx-Tb(J) 

3190 Vu( J)=f1f/<FNRhou(Tx)*Ai) 

3200! 

3210! CALCULATE INSIDE AND OUTSIDE COEFFICIENTS 
3220! 

3230 Reu)=FMRhou(Tx>=Vu)(J)«Di/FNMuw(Tx) 

3240 Cf = 1 . 

3250 0-Mf=FNCpw(Tx)*(To(J)-Ti(J)) 

3260 Qp=0/(PI*Do»L) 

3270 IF (Mi=l OR (Mi=2 AND Nrun<6)) AND J=0 THEN 
3280 Tfiln=Tsat 
3290 Kf=FMKu(Tf ilm) 

3300 Rhof =FNRhou( Tf 1 lm> 

3310 Hfg=FNHfg(Tsat)»1000 
3320 Muf=FNMuu(Tf 1 In) 

3330 Hnu=.651-(Kf '3»Rhof '2*Hfg*9.8I /(Muf*Do»0p> ) ' .3333 

3340 Tf i lnc=Tsat-0p/Hnu= .5 

3350 IF ABS((Tf ilmc-Tfilm)/Tfilmc)>.01 THEN 

3360 Tfilm=Tfilmc 



142 



3370 GOTO 3290 
3380 END IF 

3330 PRINT USING "10X.""Nasselt coefficient for first tube = "",5D.D,''“ (W/in'2> 
""" :Hnu 
3400 END IF 
3410 IF J=0 THEN 

3420 IF Mi=1 AND Nrun=1 THEN Hol=0. 

3430 PRINT " Tube Uw Heat flux Cond coef R1 R2 

RR" 

3440 PRINT •• a (tn/S) (N/m'2-K)" 

3450 END IF 
3460 f1uu=FNMuw< Tx > 

3470 Hi=FNKu(Tx)/Di'*Ci*Reu'Ex*(FNPru(Tx)) '.3333'Cf 

3480 Dt=Op/Hi*Do/Di 

3490 Cf c= < fiuw/ ( FNMuu ( Tx+D t ) ) ) ^ . 1 4 

3500 IF ftBS((Cf-Cfc)/Cfc)>.01 THEM 

3510 Cf*(Cf+Cfc)«.5 

3520 GUiU 34 /u 

3530 END IF 

3540 Lmtd=<To( J)-Ti(J))/L0G< <Tsat-Ti(J))/fTsat-To<J))) 

3550 Uo=0p/Lntd 

3560 Ho< J)=! ./< ! . /Uo-Do/(Di 'H i ) -Ru) 

3570 Rr=Uo/Ho<J> 

3580 SI=SI-^Ho(J) 

3590 IF Nrun MOD 5=1 THEN 

3600 IF Mi=1 OR <Mi=2 AND Nrun=31) THEN Ja=0 
3610 IF Mi =2 AND 5<Nrun AND Nrun<30 THEM Ja=Nrun-l 
3620 IF Mi=2 AND 35<Nrun THEN Ja=Nrun-t 
3630 END IF 

3640 IF Mi=I OR (30<Nrun AND Nrun<36 AND Mi=l) OR Nrun<6 THEN 
3650 R1=Ho<J)/Ho(0) 

3660 R2=SI/<(J+1+Ja)*Ho(0)) 

3670 ELSE 

3680 RI=Ho(J)/Ho1 

3690 R2=S1 /( ( J+1 +Ja)*Ho1 ) 

3700 END IF 
3710! 

3720! PRINT RESULTS 
3730! 

3740 PRINT USING "1 1 X . DD. 4X .DD ,DD ,2X .2(D .5DE ,2X) , 3<Z . 40 .2X) " : J+1 + Ja. J) ,0p .Ho 

(J) ,R1 .R2.Rr 

3750! 

3760! 

3770! 

3780 IF Mi=2 AND Nrun<6 AND J=0 THEN 

3790 Ho1=Ho1+Ho<0)/5 

3800 END IF 

3810 FOR K=0 TO 4 

3820 IF K=J THEN S3(K)=S3<K)+R1 

3830 IF K=J THEN S4(K) = S4(iO+R2 

3840 NEXT K 

3850 NEXT J 

3860 IF Nrun MOD 5=0 THEN 
3870 FOR K=0 TO 4 
3880 R3<K)=S3(K)/5 

3890 R4<K>=S4(K)/5 

3900 S3<K)=0. 

3910 S4(K)-0. 

3920 NEXT K 

3930 IF Mi=2 AND Nrun MOD 5=0 AND Nrun<>30 THEN As 1 =Nrun«R4(4) *Ho I 
3940! 
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oaao: rKiNi HVtrtHbc khFiU6 
3960! 

3970 PRINT " •• 

3980 PRINT " Tube 4 R3 R4" 

3990 FOR J=l TO 5 

4000 PRINT USING " 1 2X .DD .2<4X . Z . 4D) : J + Ja. R3< J- 1 ) . R4< J- 1 ) 

4010 OUTPUT i?File3:J+Ja.R3<J-1 ) ,R4< J-1 ) 

4020 NEXT J 
4030 PRINT ’• '• 

4040 END IF 

4050 IF Nrun MOD 5 = 0 AND Mi=2 AND NrunO30 AND Im=1 THEN 
4060 BEEP 

4070 PRINT USING "10X.””Set porous-tube flowmeter reading to ”",3D.D."" PERCENT 

: Set 

4080 END IF 
4090! 

4100! 

4110! 

4120 IF Im=1 THEN 
4130 BEEP 

4140 INPUT "DO YOU HAVE ANY MORE DATA ( 1 = YES . 0=NO)2” . Go_on 
4150 IF Go_on=1 THEN Repeat 
4160 ELSE 

4170 IF Mi=2 AND Nrun<30 THEN Repeat 
4180 IF Mi=1 AND Nrun<10 THEN Repeat 
4190 END IF 

4200 IF Im= 1 THEN PRINT USING ‘’lOX.DD,’"’ Data runs were stored in file , 1 0A” : 

Nrun , NewdataS 

4210 ASSIGN SFilel TO • 

4220 ASSIGN §File2 TO » 

4230 ASSIGN «File3 TO • 

4240 END 
4250! 

4260! THIS SUROUTINE CONVERTES THERMOCOUPLE VOLTAGE INTO TEMPERATURE 
4270! 

4280 SUB Tvsv(T) 

4290 COM /Cl/ C(7) 

4300 Sum=0 . 

4310 FOR 1=0 TO 7 

4320 Suin=Sujn+C( I )«T‘ I 

4330 NEXT I 

4340 T=Sum 

4350 SUBEND 

4360! 

4370! THIS FUNCTION CALCULATES PRANDTL NUMBER OF WATER IN THE 

4380! RANGE 15 TO 45 DEG C 

4390! 

4400 DEF FNPriiiYT) 

4410 Y=1 0 ' ( 1 . 03976605-T*< 1 . 3749326E-2-T*<3 . 968875E-5-3. 45026E-7*T> ) ) 

4420 RETURN Y 
4430 FNEND 

4440! 

4450! THIS FUNCTION CALCULATES THERMAL CONDUCTIVITY OF HATER 

4460! IN THE RANGE OF 15 TO 105 DEG C 

4470! 

4480 DEF FNKui(T) 

4490 Y=.5625894+T*(2.2964546E-3-T*(1 .509766E-5-4.0581652E-3'*T) ) 

4500 RETURN Y 
4510 FNEND 

4520! 
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4530! 

4540! 

4550! 

45S0 

4570 

4580 

4590 

4600! 

4610! 

4620! 

4630! 

4640 

4650 

4660 

4670 

4680! 

n o . 

4700! 

4710 

4720 

4730 

4740 

05 

4750 

4760 

4770 

4780 

4790 

4800 

4810 

4820 

4830 

4840 

4850! 

4860! 

4870! 

4880 

4890 

4900 

4910 

4920! 

4930! 

4940! 

4950! 

4960 

4970 

4980 

4990 

5000! 

5010! 

5020! 

5030! 

5040 

5050 

5060 

5070 

5080 

5090 

5100 

5110 



THIS FUNCTION CfiLCULATES SPECIFIC HEAT OF WATER 
IN THE RANGE 15 TO 45 DEG C 

DEF FNCpu<Tl 

Y=<4,21 !20358-T*<2.26826E-3-T»(4.4236lE-5+2.71428E-7*T) >)«I000 

RETURN Y 

FNEND 

THIS FUNCTION CALCULATES DENSITY OF WATER IN THE 
RANGE 15 TO 105 DEG C 

DEF FNRhow<T) 

Ro-999.52946+T*( .01269-T»<5.4825!3E-3-T*1 .234147E-5)) 

RETURN Ro 
FNEND 



TtlT^ r-Ifkir'TTOM 

1(1X0 I I XOM 



APPLIES 



CORRECT 



TOkK"' 
X oi «o 



TO 



THERMOCOUPLE READINGS 



DEF FNTemp<T.I) 

DIM A(14) .8(14) 

DATA 0 . 640533 . 0 .573054 .0.593101.0. 57298 . 0 . 56228 . 0 . 567384 . 0 . 569577 

DATA 0.553951 .0.552008.0.566955,0.520998.0.522661 .0.531008,0.560788.0.5524 

DATA 1 1 .8744.8.63163.9.39412,3.570246.3.299436.8.36677.8.04507.7.459766 
DATA 7 . 498928 . 7 . 9408 , 5 . 87072 , 5 . 39 1 556 , 6 . 1 3399 . 6 . 48586 . 6 . 326224 
READ A(*) ,B<«> 

IF K15 THEN 
T=T-<A<I)-B<I)*.001»T) 

ELSE 
T=T-.5 
END IF 
RETURN T 
FNEND 

THIS FUNCTION COMPUTES THE SPECIFIC VOLUME OF STEAM 



DEF FNVvst(T) 

V=58.4525588-T*< 1 .5 I508776-T'*< .01372746585-T-4 .25366711E-5) > 

RETURN V 

FNEND 



THIS FUNCTION CONVERTS THE VOLTAGE READING OF THE PRESSURE 
TRANSIDUCER INTO PRESSURE IN MM HG 



DEF FNPvsv<V> 

Y=1 .1 1 03462+ 163. 3641 3*V 

RETURN Y 

FNEND 



THIS FUNCTION CALCULATES THE SATURATION TEMPERATURE OF STEAM AS A FUNCTION 
OF PRESSURE 



DEF FNTvsp<P) 

IF P<600 THEN 

T=31 .8776153+P*< .235854929-P*(3.6613664E-(i-P»2.4I652372E-7)) 
ELSE 

T-59.3G562+P*( .07379467-P*<3. 15662E-5-P«6.27246E-9>) 

END IF 
RETURN T 
FNEND 
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5!20! 

5130! THIS FUNCTIOM COMPUTES THE VISCOSITY OF WATER 
5140! 

5150 DEF FNMuui(T> 

5160 Mu»! .57609473E-3-T»(3.51 198576E-5-T*<3.5835816E-7-1 .3655861 15E-9*T)> 
5170 • RETURN Ma 
5180 FNEND 
5190! 

5200! THIS FUNCTION COMPUTES THE LATENT HEAT OF VAPORIZATION 
5210! 

5220 DEF FNHfg(T> 

5230 Hfg=2497.7389-T*(2.2074+T*<1 .7079E-3-2.8593E-6»T>) 

5240 RETURN Hfg 
5250 FNEND 
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MODIFIED WILSON-PLOT METHOD 



1 000 ! 

lOloi FILE NAME: WILSON 
1020 ’ 

1030! THIS PROGRAM COMPUTES THE SIEDER-TATE 
1040! COEFFICIENT FOR FLOW IN TUBES 
1 050 ! **•*»**»»»*****»* 

1060 COM /Cc/ C<7) 

1070 DIM Vui( 18) .Tci< 18) ,Tco( 18) .l3< 18) .Md< 18) .Ta(18) 

1080 DATA 0.10085091 .25727.94369.-767345.8295.78025595.81 
1090 DATA -9247486589. 6. 97688E+1 1 , -2. 6G192E+13. 3. 94078E+I4 
1100 READ C(») 

1110 PRINTER IS 701 
1120 BEE^ 

1130 CLEAR 709 

1140 INPUT "ENTER MONTH. DATE AMD TIME (MM : DD :HH:MM :SS)" . BS 
1150 OUTPUT 709:"TD":BS 
1160 Di=.0141 

1170 Ai=PI*^Di'2/^ 

1130 Do=. 015875 
1190 L=.305 

1200 Km=21.9 

1210 Ru=(Do-Di )»Do/(Km'»(Do + Di ) ) 

1220 Rfi=0. 

1230 Senes: ! 

1240 OUTPUT 709;"TD" 

1250 ENTER 709; AS 

1260 PRINT USING ”1 OX ."“Month, date .and time: ”".15A":A$ 

1270 BEEP 

1280 INPUT "ENTER INITIAL GUESS FOR SIEDER-TATE COEFFICIENT" .Ci 
1290 BEEP 

1300 INPUT "ENTER E-XPCNEMT FOR REYNOLDS NUMBER". Xn 

1310 PRINT USING "1 OX . ""Ini t lal guess for Sieder-Tate coefficient = ”".Z.DD";Ci 
1320 PRINT USING "1 OX , ""Exponent for the Reynolds n<imber = ”",2.DD";Xn 

1330 BEEP 

1340 INPUT "ENTER THE INPUT MODE < 1=3054A .2=FILE)" . Im 
1350 IF Im=1 THEN 
1360 BEEP 

1370 INPUT "GIVE A NAME FOR THE DATA FILE" . D_f i leS 
1380 CREATE BDAT D_fi!e$,5 
1390 ELSE 
1400 BEEP 

1410 INPUT "GIVE THE NAME OF THE DATA F ILE" . D_f i leS 

1420 PRINT USING "1 CX , "”Fol ioying analysis uas performed for data stored in fii 
e 10A":D_f 1 leS 
1430 BEEP 

1440 INPUT "ENTER THE NUMBER OF DATA RUMS STORED". Nrun 
1450 END IF 

1460 ASSIGN SFile TO D_file$ 

1470 BEEP 

1480 INPUT "GIVE A NAME FOR PLOTTING DATA FILE",DolotS 
1490 CREATE BDAT DplotS.5 
1500 ASSIGN JiFilep TO DplotS 
1510 BEEP 

1520 INPUT "ENTER ANALYSIS TYPE < 1 -HI . 2*UI ) " . I t 

1530 IF It-1 THEN PRINT USING "1 OX . ""Analys is type ■ Hi-METHOD 

1540 IF It=2 THEN PRINT USING "1 OX . ""Analysis type - Ui-METHOD 
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1550 K=0 

1560 PRINT •• •• 

1570 PRINT •• Data Vu Tsat Ti To’* 

1580 PRINT " f <m/S) (C) (C) (C)” 

1 590 Repeat : ! 

1600! 



1610! 

1620! 

1630! 

1640! 

1650 

1660 

1670 

1680 

1690 

1700 



RECORDS THERMOCOUPLE ftND PRESSURE TRANSDUCER 
READINGS AUTOMATICALLY THROUGH THE HP 3054A 
AUTOMATIC DATA ACQUISITION/CONTROL SYSTEM 



IF Im=l then 
BEEP 

INPUT "ENTER 



OUTPUT 709 
OUTPUT 722 
OUTPUT 709 



"AR 

”F1 

"AS 



FLOWMETER READING (AS A PERCENT )" .Fm 



AFO ALO" 
R1 TI Z1 
SA" 



FLl' 



i7;01 READS THERMOCOUPLE FOR MATPR jnlET 
1720 ENTER 722;T<0) 

1730 OUTPUT 709:"AR AF3 AL5" 

1740 OUTPUT 722; "FI R1 TI Z1 FLl" 

1750 FOR I-l TO 3 
1760 OUTPUT 709: "AS SA" 

1770! READS THREE THERMOCOUPLES FOR WATER OUTLET 
1780 ENTER 722:T<I) 

1790 NEXT I 

1800 OUTPUT 709:"AR AF19 AL19" 

1810 OUTPUT 722:"F1 R1 TI Z1 FLl" 

1820 OUTPUT 709: "AS SA” 

1830! READS PRESSURE TRANSDUCER 
1840 ENTER 722:P_volts 
1850 Pvap=FNPvsv ( P_vo 1 ts ) 

1860 Ts(K)=FNTvsp<Pvap) 

1870 Sujn-O. 

1880 FOR 1=0 TO 3 

1890! CONVERT VOLTAGE READINGS TO TEMPERATURE 
1900 CALL Tvsv<T(I)> 

1910 IF 1=0 THEN 

1920 Tci <K) =FNTefnp< T ( 0 ) . 0 ) 

1930 ELSE 
1940 M=I+2 



1950! APPLY THERMOCOUPLE CORRECTIONS 
i960 To=FNTenp< T( I) ,M> 

1970 Sujn=Sujii+To 
1980 END IF 
1990 NEXT I 
2000 Tco<K)=Siun/3. 

2010 ELSE 

2020 ENTER 3F i le : Ts(K) . Tci (K) . Tco(K) ,Fm 
2030 END IF 

2040 Ta(K)=(Tci<K)+Tco(K))*.5 

2050 Md<K)=66.86»Fm/(100»60) 

2060! COMPUTE WATER-SIDE VELOCITY 
2070 Vui ( K ) =Md ( K) / (FNRhow ( Ta( K) ) *A i ) 

2080 BEEP 

2090 IF Im=l THEN INPUT "ARE YOU TAKING MORE DATA < I =YES . O-NO)’?" . Go 
2100 M=K+1 

2110 PRINT USING "1 OX .DD . 5(2X .DDD . DO) ” :M . Vu<K) . Ts(K) . Tci <K) . Tco(K> 
2120 IF Itn=l THEN OUTPUT fiFile:Ts(K) ,Tci(K) ,Tco<K) .Fm 
2130 K=K+1 

2140 IF Im-1 THEN 



on 
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2!50 IF Go on=l THEN Repeat 
2160 ELSE 

2170 IF fKNrun THEN Repeat 

2180 END IF 

2190 K=K-I 

2200 J=0 

2210 Jj=0 

2220? 

2230? PERFORM ITERATION TO COMPUTE SIEDER-TATE COEFFICIENT 
2240? 

2250 Ssq=0 
2260 Sx=0. 

2270 Sy=0. 

2280 Sxs»0. 

2290 Sxy=0. 

2300 J=J+1 

2310 PRINT •• •• 

2320 PRINT ” Iteration number = ":J 

2330 IF J=1 OR Jj=l THEN PRINT ” X 1/HI <TW-TB) Q 

CF" 

2340 FOR 1=0 TO < 

2350 0=Md<I)*(Tco(I)-Tci<I))*4180. 

2360 Lr>td=<Tco<I)-Tci(I))/LOG(<Ts<I)-Tci(I))/(Ts<I)-Tco<I) )) 

2370 Un=Q/(Lmtd»PI*Do«L> 

2380 Unr=1./Un 

2390 Rei=FNRhoui<Ta< I) )->Vu(I)*Di/FNMu(Ta<I)> 

2400 Prw=FNPru.'( Ta< I ) ) 

2410 X = Rei ' ( -Xn)»Pru ■ (- .3333) 

2420 IF It»2 THEN 
2430 H i r=Unr»D 1 /Do 
2440 Ku=FNKu(Ta(I)) 

2450 GOTO 2660 
2460 END IF 
2470 IF 1=0 THEN 
2480 Cf=1. 

2490 Ku=FNKuj<Ta<I) ) 

2500 Hi-Kw/Di*Ci»RerXn*Pru' .3333*Cf 
2510 Dt = Q/(PI-‘Di*L»Hi) 

2520 Cfc»(FNMu(Ta<I) )/FNMa<Ta(I)+Dt))' .14 
2530 IF ABS<<Cfc-Cf)/Cfc)>.01 THEN 
2540 Cf=(Cf+Cfc)/2. 

2550 GOTO 2500 
2560 END IF 

2570 Ho-I ./(Unr-Ru-Rf i»Do/Di-Do/(Di«Hi)> 

2580 Hir=1/Hl 
2590 0o=0 

2600 ELSE 
2610? COMPUTE 1/HI 

2620 Hir=<Unr-i/Ho*<0/Qo)'. 3333*Rw*Rf i«Do/Di ) »D i /Do 

2630 Dt-0/(PI*‘Di«L>*Hir 

2640 Cf = < FNMu< Ta< I ) ) /FNMu< Ta< I ) +D t ) ) '' . 1 4 

2650 END IF 

2660 IF It=2 THEN 

2670 Cf = 1 . 

2680 Hic=Ku*Ci/Di-Rel‘Xn«Pr.j-'.3333‘Cf 
2690 Dt=0/<PI*Di»L*Hic) 

2700 Cfc=(FNMa(Ta(I))/FNMa<Ta(I)+Dt))".14 
2710 IF ABS< <Cf-Cfc)/Cfc)> .01 THEN 
2720 Cf»<Cf+Cfc)«.5 
2730 GOTO 2680 
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2740 END IF 
2750 END IF 
2760 Op'Q/(PI*Do*L) 

2770 X=X/Cf 

2780 IF Jj=1 THEN 

2790 OUTPUT @Filep;X,Hir 

2800 Hirc=A+B«X 

2810 Ssq-Ssq+<Hir-Hirc)"2 

2820 END IF 

2830 IF J=1 OR Jj=t THEN PRINT USING " 1 0X . Z . 7D . 3X . Z . 7D . 4X ,DD . DD .2X . D . 3DE . X . Z . 5D 
•':X,Hir ,Dt .Op.Cf 

2840! COMPUTE COEFFICIENTS FOR LEAST-SQUARES SCHEME 

2850 Sx=Sx+X 

2860 Sy=Sy+Hir 

2870 Sxs-Sxs+X'2 

2880 Sxy=Sxy+X*Hir 

2890 NEXT T 

2900 N=K+1 

2910! COMPUTE THE SLOPE OF THE LEAST-SQUARES LINE 
2920! DEVELOPED FOR THE WILSON PLOT 
2930 B= (N»Sxy-Sy»Sx ) / < N*Sxs-Sx*Sx ) 

2940 Cic=Di/(B*Ku) 

2950 PRINT USING "I OX . ""Intermediate value of Sieder-Tate coefficient = "",Z.4D 
":Cic 

2960 IF ABS((Ci-Cic)/Cic)>.01 THEN 

2970 Ci=Cic 

2980 GOTO 2260 

2990 END IF 

3000 IF Jj=! THEN 

3010 PRINT USING " 1 OX . ""Sieder-Tate coefficient = "",Z.4D";Cic 

3020 A=(Sy-B»Sx)-/M 

3030 PRINT USING "1 OX . ""Estimated fouling factor = "".MZ.5DE."" (m'2-K/W) :A 

3040 PRINT " Least-squares line:" 

3050 PRINT USING "13X. ""Slope = "".Z.5D";B 

3060 PRINT USING "13X . ""Intercept = "",MZ.5DE":A 

3070 ELSE 

3080 Jj=1 

3090 GOTO 2260 

3100 END IF 

3110 PRINT USING "i OX. ""Sum of squares = . D . 5DE" :Ssq 

3120 ASSIGN '"'File TO • 

3130 ASSIGN «Fiiep TO » 

3140 IF Im=l THEN 
3150 BEEP 

3160 PRINT USING "10X.""NOTE: "",DD."" Data runs were stored in file "".8A":M,D 
_f i 1 e$ 

3170 END IF 

3180 PRINT USING "10X.""N0TE: "".DD."" X-Y pairs were stored in file "".10A":Ki- 
1 . Dp lots 
3190 BEEP 

3200 INPUT "ARE YOU RUNNING ANOTHER SERIES < 1 ' YES . 0=NO>2" .Go_on 
3210 IF Go_on=1 THEN Series 
3220 END 

3230 DEF FMRhow(T) 

3240! 

3250! THIS FUNCTION COMPUTES THE DENSITY OF WATER 
3260! 

3270 Ro=1006.35724-T-< .774489-T-(2.262459E-2-T*3.0330qE-4>) 

3280 RETURN Ro 
3290 FNEND 
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3300 DEF FNPvsv<V) 

3310! 

3320! THIS FUNCTION CONVERTS THE PRESSURE TRANSDUCER 
3330! READING FROM VOLTS TO PRESSURE IN MM HG 
3340! 

3350 Y-1 . 1103462+163. 364I3*V 

3360 RETURN Y 

3370 FNEND 

3380 DEF FNTvsp(P) 

3390!- 

3400! THIS FUNCTION COMPUTES THE SATURATION TEMPERATURE 

3410! CORRESPONDING TO PRESSURE IN MM HG 

3420! 

3430 IF P<600 THEN 

3440 T=31 .8776158+P*< .235854929-P*(3.6613664E-4-P*2.41652372E-7) ) 

3450 ELSE 

3460 T=59.3G562+P»( .073794b/-P*(3. I5662E-5-F-S.27246E-9) ) 

3470 END IF 
3480 RETURN T 
3490 FNEND 
3500 DEF FNTet7>p<T.I) 

3510! 

3520! THIS FUNCTION APPLIES THERMOCOUPLE CORRECTIONS 
3530! 

3540 DIMA(5).B(5) 

3550 DATA 0 .640533 . 0 . 573054 . 0 .5931 0 1 . 0 . 57298 . 0 . 567384 .0 .569577 
3560 DATA 1 1 . 8744 . 8 . 631 63 .9 . 3941 2 .8 .570246 , 8 . 299436 .8 . 36677 
3570 READ AC) .B(*> 

3580 T = T-(A<I)-B< D-.OOCT) 

3590 RETURN T 
3600 FNEND 
3610 SUB Tvsv(T) 

3620 COM /Cc/ C<7) 

3630 Sujn=0 . 

3640 FOR 1=0 TO 7 

3650 Sam=Sun+C< I ) '•T ' I 

3660 NEXT I 

3670 T=S.im 

3680 SUBEND 

3690! 

3700! THIS FUNCTION COMPUTES PRANDTL NUMBER FOR WATER 
3710! 

3720 DEF FNPr-j<T) 

3730 Pr=10 •< : .09976605-T-< .013749326-T+<3.968875E-5+3.45026E-7*T) ) ) 

3740 RETURN Pr 
3750 FNEND 
3760 DEF FNMa<T) 

3770! 

3780! THIS FUNCTION COMPUTES THE VISCOSITY OF WATER 
3790! 

3800 Mu=1 .5037546575E-3-T*^<3.025732489E-5-T><2.626439826E-7-T-8.18601937E-10)) 

3810 RETURN Mu 

3820 FNEND 

3830 DEF FNKu(T) 

3840! 

3850! THIS FUNCTION COMPUTES THE THERMAL CONDUCTIVITY OF WATER 
3860! 

3870 Kuia= .572183504477 + 1 .52770121209E-3-T 
•3880 RETURN Kwa 
3890 FNEND 
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C. PLOTTING PROGRAM 



1000! FILE NAME: PLOT 
1010 PRINTER IS 705 
1020 BEEP 

1030 INPUT “ENTER MINIMUM AMD MAXIMUM X-UALUES” .Xnin.Xmax 
1040 BEEP 

1050 INPUT “ENTER MINIMUM AMD MAXIMUM Y-UALUES“ . Ynin. Ymax 
1060 BEEP 

1070 INPUT “ENTER STEP SIZE FOR X-AXIS" . Xstep 
1080 BEEP 

1090 INPUT "ENTER STEP SIZE FUR Y-AXIS" . Ystep 
1100 BEEP 

1110 PRINT “IN:SP1;IP 2300.1800.8300.6800:" 

1120 PRINT "SC 0.100.0.100:TL 2.0:" 

Ilou OIA-IUU/' AriCl-A “An X ll / 

1140 Sfy=1 00/< Ymax-Ynin) 

1150 PRINT "PU 0.0 PD" 

1160 FOR Xa=Xrun TO Xmax STEP Xstep 

1170 X=(Xa-Xmin)*Sfx 

1180 PRINT "PA":X.".0: XT:" 

1190 NEXT Xa 

1200 PRINT "PA 100.0:PU;“ 

1210 PRINT "PU PA 0.0 PO" 



1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 



FOR Ya=Ymin TO Ymax STEP Ystep 

Y-(Ya-Ymin)*Sfy 

PRINT "PA 0.":Y."YT" 

NEXT Ya 

PRINT "PA 0.100 TL 0 2" 

FOR Xa=Xnin TO Xmax STEP Xstep 

X=(Xa-Xnin)«Sfx 

PRINT "PA":X,".100; XT" 

NEXT Xa 

PRINT "PA 100.100 PU PA 100.0 PD" 

FOR Ya=Ymin TO Ymax STEP Ystep 

Y=(Ya-Ymin)*Sfy 

PRINT "PD PA 100.".Y."YT" 

NEXT Ya 

PRINT "PA 100.100 PU" 

PRINT "PA 0.-2 SR 1 .5.2“ 

FOR Xa=Xnin TO Xmax STEP Xstep 
X=<Xa-Xnin)»Sf x 
PRINT “PA“:X.".0:“ 

PRINT "CP -2.-1 :LB":Xa:““ 

NEXT Xa 

PRINT "PU PA 0.0" 

FOR Ya-Ymin TO Ymax STEP Ystep 

Y-(Ya-Ymin)*Sfy 

PRINT "PA 0.":Y.“" 

PRINT "CP -4,-.25:LB":Ya:"“ 

NEXT Ya 
BEEP 

INPUT "ENTER X-LABEL" .XlabelS 
BEEP 

INPUT "ENTER Y-LABEL" . YlabelS 

PRINT "SR I.5.2:PU PA 50.-10 CP" : -LEM( XlabelS) /2: "0 :LB" 
PRINT "PA -11.50 CP 0.";-LEN<YlabelS)/2*5/6:"DI 0.1 :LB" 
PRINT "CP 0.0 DI" 



XlabelS 
Y labels 
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1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

16<!i0 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

173C 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2030 

2090 

2100 

2110 

2120 

2130 

2140 

2150 



Repeat : ! 

BEEP 

INPUT "DO YOU W8MT TO PLOT DATA FROM A FILE ( t=YES.O=NO)?".Ok 

IF Ok=1 THEN 

BEEP 



INPUT "ENTER THE NAME OF THE DATA FILE" .D_f i leS 

ASSIGN ^File TO D_fileS 

BEEP 

INPUT "ENTER THE BEGINNING RUN NUMBER". Md 
BEEP 

INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED" .Npai rs 
BEEP 

INPUT "SELECT A SYMBOL FOR THE PLOTTER ( 1 =* .2=+ .3=c . 4=o .5= 
PRINT "PU DI" 

IF Sy=l THEN PRINT "SM»" 

IF Sv»2 THEN PRINT "SM+" 

IF Sy-3 THEM PRINT "SMc" 

Tfr r* _ /. Trtrii PDTMT 
ir rr\li«( 

IF Sy=5 THEM PRINT "SM'" 



BEEP 

INPUT "SELECT MODE 
IF Md>Noairs THEN 
FOR 1=1 TO (Md-I ) 
IF Jn=1 THEN ENTER 
IF Jn=2 THEN ENTER 
NEXT I 
END IF 

FOR 1=1 TO Npairs 
IF Jm=! THEN ENTER 
IF Jni = 2 THEN ENTER 



( I’HN/HI .2=HN(avg)/H1 )". Jm 



MiFi ie :Xa. Ya, Yy 
i<File:Xa. Yy ,Ya 



^iF 1 le : Xa, Ya. Yy 
!sFile:Xa.Yy.Ya 



X=<Xa-Xnin)-*Sf X 



)".Sy 



Y=(Ya-Ymin>»Sfy 
PRINT "PA".X.Y."" 
NEXT I 
BEEP 



ASSIGN 5Flle TO * 

INPUT "DO YOU HAVE MORE DATA TO BE PLOTTED (I=YES.O=NO)?" .Go_on 
IF Go_on=1 THEN Repeat 
END IF 
BEEP 

INPUT "DO YOU LIKE TO PLOT THE NUSSELT RELATION ( 1 =YES . 0=NO)2" . Go_on 
PRINT "PU:SM" 

IF Go_on=l THEN 

FOR Xa=Xmin TO Xmax STEP Xstep/50 
X=< Xa-Xmn)*'Sf x 

IF Jni=l AND Xa>Xmin THEN Ya=Xa' . 75-< Xa- 1 ) ' . 75 
IF Jm=2 AND Xa>Xirun THEN Ya=Xa'(-.25> 

IF Xa=Xmin THEN Ya=1 
Y=(Ya-Ymin)»Sfy 
PRINT "PA".X.Y,"PD" 

NEXT Xa 
BEEP 

PRINT "PU" 

INPUT "MOVE THE PEN TO LABEL THE NUSSELT LINE". Ok 



PRINT "LBNusselt" 

END IF 
BEEP 

INPUT "DO YOU LIKE TO PLOT KERN RELATIOMSH IP?" . Yes 
IF Yes=1 THEN 

tOR Xa=Xmn TO Xmax STEP Xstep/20 
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21S0 Ya-Xa (-1/6) 

2170 X=(Xa-Xmin)^Sfx 
2180 Y=(Ya-YiTun)*Sfy 
2190 PRINT ••PA”.X. Y."PD'’ 

2200 NEXT Xa 
2210 PRINT "PU" 

2220 ■ BEEP 

2230 INPUT “MOV/E THE PEN TO LABEL KERN RELATIONSHIP" .Ok 
2240 PRINT "LBKern:PU" 

2250 END IF 

2260 PRINT "PU PA 0.0" 

2270 BEEP 

2230 INPUT "OK TO PLOT CURVE FOR NEW DESIGN ( 1 =0K .0=N0)7" .Ok 
2290 IF Ok=1 THEN 

2300 FOR Xa=Xnin TO Xmax STEP Xstep/2 
2310 IF Xa=1 THEN Ya=l 

2320 IF Xa>i THEN Ya= < < Xa»5) ' . 75-<Xa«5- 1 ) ' . 75+5) /6 
2331) X=(Xa-Xmin)»Sf X 
2340 Y=(Ya-Yrun)->Sfy 
2350 PRINT "PA".X. Y."PD" 

2360 NEXT Xa 
2370 PRINT "PU" 

2380 END IF 
2390 BEEP 

2400 INPUT "DO YOU LIKE TO PLOT EISSENBERG RELATION^" ,Go_on 

2410 IF Go on=1 THEN 

2420 FOR Xa=Xiran TO Xmax STEP Xstep/10 

2430 Ya-.6+.42»Xa'-(-.25) 

2440 X=<Xa-Xmin)*Sfx 
2450 Y=(Ya-Ymin)-Sfy 
2460 PRINT ••PA".X.Y.."PD" 

2470 NEXT Xa 
2480 PRINT "PU" 

2490 BEEP 

2500 INPUT "MOVE THE PEN TO LABEL THE EISSENBERG LINE". Ok 
2510 PRINT "LBEissenberg:PU" 

2520 END IF 
2530 PRINT "PU" 

2540 BEEP 

2550 INPUT "DO YOU LIKE TO PLOT THE EXPTL CURVE ( 1 =Y . 0=N0)?" . Go_on 
2560 IF Go_on=1 THEN 
2570 BEEP 

2580 INPUT "ENTER THE EXPONENT". Ex 
2590 FOR Xa=Xmin TO Xmax STEP Xstep/10 
2600 Ya=Xa'(-Ex) 

2610 X=(Xa-Xmin)*Sfx 
2620 Y=(Ya-Ynin)*Sfy 
2630 PRINT "PA".X. Y."PD". 

2640 NEXT Xa 
2650 PRINT "PU" 

2660 BEEP 

2670 INPUT "MOVE THE PEN TO LABEL THE EXPTL CURVE". Ok 
2680 PRINT "LBHN<avg>/H1 =M" 

2690 PRINT "PR 1.1" 

2700 PRINT "LB": -Ex:"" 

2710 END IF 
2720 BEEP 

2730 INPUT "DO YOU LIKE TO DRAW A STRAIGHT LINE7".Go_on 
2740 IF Go_on-1 THEN 
2750 BEEP 

2760 INPUT "ENTER THE SLOPE", SI 
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2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2390 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 



BEEP 

INPUT "ENTER THE INTERCEPT" .He 

FOR Xa=Xmin TO Xmax STEP <Xmax-Xmin> 

Ya=flc+Si»Xa 

Y’<Ya-Ymin)»Sfy 

X = <Xa-Xrun)’»Sfx 

IF Y<0 THEN 

Xatn» f Ym i n-Ac > /S 1 

X=<Xatn-Xinin)'»Sfx 

Y=0 

END IF 

IF Y>100 THEN 
Xam= ( Ytnax-ftc / /S 1 
X = ( Xar>-Xnun) «Sf x 
Y = 100 
END IF 

PRINT "PA".X. Y."PD" 

NEXT Xa 
END IF 

PRINT "PU SPO" 

END 



155 



LIST OF REFERENCES 



1. Standards for Steam Surface Condensers ^ 6th ed.. 

Heat Exchange Institute, 1970. 

2 . Standards of_Tubular_ Exchange Manufacturers 

Association , 4th ed.. Tubular Exchanger Manufacturers 
Association. Inc., 1959. 

3. Search, H.T., A Feasibility Study of Heat Tra nsfer 
Improvement, in Marine Steam Condensers . MSME, Naval 
Postgraduate School, Monterey. CA, January 1977. 

4. Webb, R.L., "The Use of Enhanced Surface Geometries in 

Condensers: An Overview," Power Condenser Heat 

Transfer Technology . Marto, P.J. and Nunn, R.H. 

(Eds), Hemisphere Publishing Corp., New York, 1980, 
pp. 287-318. 

5. Beck, A.C., A Test Facility to Measure Heat Transfer 
Performance of Advanced Condenser Tubes . MSME, Naval 
Postgraduate School. Monterey, CA, March 1978. 

6. Pence, D.T., An E xp erimental Study of Steam Condensa - 
tion on a Single Horizontal T u be , MSME, Naval 
Postgraduate School, Monterey, CA, March 1978. 

7. Reilly, D.J. , An Experimental investigation. on 
Enhanced. Heat Transfer on J)PXi.zpri:taI_C c>ndeh5 g.i:_ T u b g$, 
MSME, Naval Postgraduate School. Monterey, CA, 

March 1978. 

8. Fenner, J.H., An Exp er i mental C omparison of Enhanced 
Heat Transfer Condenser Tub ing, MSME. Naval Post- 
graduate School, Monterey, CA. September 1978. 

9. Ciftici, H. , An Experimental Study of Filmwise 
Condensation, on_ Horizontal Enhanced Condenser Tubing , 
MSME, Naval Postgraduate School. Monterey, CA, 

December 1979. 

10. Marto, P.J. , Reilly, D.J. , Fenner, J.H., "An 
Experimental Comparison of Enhanced Heat Transfer 
Tubing," Advances in Enhanced _Heat Transfer . 

Proceeding of the 18th National Conference. American 
Society of Mechanical Engineers. New York, NY, 1979. 



156 



11. 



Eissenberg, D.M., An In vestigation of the Variables 
Aff.ectinq_ Steam Condensation on the Outside .of _a 
Horizontal Tube Bundle . Ph.D Thesis- University of 
Tennessee. Knoxville- TN, December 1972. 

12. Noftz, P.J. , Effects of Condensate Inundation and 
Vapor Velocity on Heat Transfer in a Condenser Tube 
Bundle . MSME, Naval Postgraduate School. Monterey. 

CA. June 19^82. 

13. Morrison, R.H., A Test Condenser to Measure Condensate 
Inu ndation Effects in a Tube Bundle . MSME. Naval 
Postgraduate School. Monterey, CA, March 1981. 

14. Thomas, A., Lorenz, J. J. , Hillis- D.A. , Young, D.T. and 
Sather. N.F., "Performance Tests of the 1 Mwt Shell and 
Tube Exchangers for OTEC", Proc. 6th OTEC Conf. . 

Paper 3c. June 19-22, 1979. 

15. Cunningham, J. , "The Effect of Noncondensable Gases on 
Enhanced Surface Condensation." Power Condenser Heat 
Transfer Technology . Marto, P.J. and Nunn, R.H. 

(Eds), Hemisphere Publishing Corp. New York, 1980, 

pp. 353-365. 

16. Newson, I.H., "Enhanced Heat Transfer Condenser Tubing 
for Advanced Multistage Flash Distillation Plants", 

5th Int. Symp. on Fresh Water from Sea . Vol. 2, 

pp. 107-115, 1976. 

17. Eissenberg, D. and Bugue . Do. "Tests on an Enhanced 
Horizontal Tube Condenser Under Conditions of 
Horizontal Steam Cross Flow," 4th Int. Heat Transfer 
Conf. , Versailles- H.E. 2, I, September 1976. 

18. Nobbs, D.W. , The Effect of Downward Vapour Velocity 
and Inundation on the Condensate Rates on Ho rizontal 
Tubes and Tube Banks . Ph.D Thesis, University of 
Bristol- Bristol, England. April 1975. 

19. Jakob, Max, H eat Transfer . Vol. 1, J. Wiley & Sons, 
Inc., 1948, pp. 667-673. 

20. Kern, D.Q., "Mathematical Development of Loading in 
Horizontal Condensers." AICHE Journal . Vol. 4, 1958, 
pp. 157-160. 

21. Chen. M.M., "An Analytical Study of Laminar Film 
Condensation; Part 2 - Single and Multiple Horizontal 
Tubes." Journal of He at Tr ansfer . February 1965. 



157 



22, Berman. L,D. , "Heat Transfer with Steam Condensation on 
a Bundle of Horizontal Tubes,” Thermal Engineering , 

Vol. 28, 1981, pp. 218-224. 

23. Naval Postgraduate School Report NPS69-82-005, 

Computer Rating of Enhanced .Marine Condensers , by 
R.H. Nunn and P.J. Marto, p. 13, August 1982, 

24, Berman, L.D, and Tumanov, Y.A. , "Investigation of Heat 
Transfer in the Condensation of Moving Steam on a 
Horizontal Tube," Teoloenergitika , Vol. 9, 1962, 

pp. 77-83. 

25, Eissenberg, D.M., "Combined Effects of Vapor Shear and 
Inundation in Tube Bundles." Power Condenser Heat 
Transfer Technology . Marto. P.J. and Nunn, R.H, 

(Eds) , Hemisphere Publishing Corp, New York, 1980, 
pp. 225-227. 

26. Metals Handbook . 8th Ed., Vol 2, p. 665, American 
Society for Metals, 1964. 

27, Holman. J.P., H eat Transfer . 4th Ed., McGraw-Hill, 

1976. 

28. Wilson, E.E., A Basis fo r Rational Design .of _Heat 
Transfer Apparatus , paper presented at the Spring 
Meeting of the Society of Mechanical Engineers. 

Buffalo. NY, June 1935. 

29, Wanniarachchi , S.A. , "Modified Wilson Plot Technique," 
private communication. January-June 1983. 

30. Incropera, P.F. and DeWitt. P.D. , Fund amentals of Heat 
Transfer . John Wiley & Sons. 

31. Kline. S.J, and McClintock. F.A. , Describing 
U ncertainties in Single Experiments . Mech. Eng in., 

Vol, 74, pp. 3-8, January 1953, 



158 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Technical Information Center 2 

Cameron Station 

Alexandria, Virqinia 22314 

2. Library, Code 0142 2 

Naval Postgraduate School 

Monterey. California 93940 

3. Department Chairman. Code 69 1 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, California 93940 

4. Professor P.J. Marto, Code 69Mx 5 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey. California 93940 

5. Professor R.H. Nunn, Code 69Nn 1 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey. California 93940 

6. Dr- A.S- Wanniarachchi , Code 69 1 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, California 93940 

7. Mr. R-W. Kornbau 2 

Code 2721 



David Taylor Naval Ship Research and 
Development Center 
Bethesda. Maryland 20084 

8. Hellenic Navy General Staff 3 

Department of Education 
% Embassy of Greece 
Office of Naval Attache 
2228 Massachusetts Avenue. N.W. 

Washington, DC 20008 



159 



9. Professor Th. Boufounos 1 

Naval Academy (S.N.D.) 

% Embassy of Greece 

Office of Naval Attache 

2228 Massachusetts Avenue, N.W. 

Washington, DC 20008 

10. Professor G. Malachias 
Naval Academy (S.N.D.) 

% Embassy of Greece 
Office of Naval Attache 
2228 Massachusetts Avenue, 

Washington, DC 20008 

11. Professor I. loannides 1 

Department of Naval Engineering 

National Technical University of Athens 
42, October 28th Street 
Athens 147 . Greece 

12. LT Georgies D, Kanakis, H.N. 3 

2 Ersis, Paleon Psychicon 

Athens, Greece 

13. LT loannis Georgiadis, H.N. 1 

SMC 3001 

Naval Postgraduate School 
Monterey. California 93940 



1 



N.W. 



160 



■i 








201G88 



Thesis 

K1325 Kanakis , 

c.l The effect of conden-^^ 

sate inundation on 
steam condensation heat 
transfer to wire- 
wrapped tubing. 

1 7 



n 

a heat 



26 mr 



3 7 7 1 7 



201638 

Thesis 

K1325 Kanakis 

c.l The effect of conden- 

sate inundation on 
steam condensation heat 
transfer to wire- 
wrapped tubing. 



